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ABSTRACT 

Useful  estimates  of  seabed  acoustic  impedance  variations 
in  depth  and  location  can  be  derived  from  measurements  of  the 
travel  time  rate  at  which  low  frequency  normal  incidence  re- 
flected seismic  energy  is  received.  Analysis  techniques  are 
straightforward  and  easily  applied  to  station  data  and  to  tra- 
verse profiles  to  evaluate  the  lateral  continuity  of  material 
properties  from  one  location  to  another.  Added  to  critical 
judgment,  geological  and  geophysical  evidence,  e.g.  drill-hole 
cores,  outcrops  and  interval  and  critical  refraction  velocity 
measurements,  reflection  energetics  measurements  provide  im- 
proved insight  with  which  to  justify  classifications  of  seabed 
materials  within  the  constraints  of  available  information. 
Methods  to  map  seabed  response  characteristics  are  discussed. 

The  analysis  technique  is  tested  using  reflection  data 
from  six  stations  and  two  profiles  in  the  environs  of  the  Somali 
Basin.  Two  clearly  different  impedance-depth  distributions  are 
evident,  one  to  the  northwest  of  Chain  Ridge  is  near- linear  down 
through  the  sediments  and  including  basement,  the  other  to  the 
south  of  the  basin  proper  exhibits  a near-uniform  low  impedance 
throughout  the  sediment  column  with  sharp  increases  at  or  near 
basement,  in  agreement  with  that  found  nearby  at  Deep  Sea  Dril- 
ling Project  site  #235.  Continuity  of  material  properties  is 
more  uniform  in  the  southern  data  than  the  northern  where  local 
variations  in  seafloor  reflectivity  and  seabed  response  are  evi- 
dent. 


Major  conclusions  from  our  "audit"  of  reflected  energy  are 
that  reflection  energetics  measurements  and  impedance  estimates 
agree  with  and/or  corroborate  the  other  geological  and  geophysical 
evidence,  and  that  such  an  audit  has  unique  utility  in  the  inter- 
pretation of  seismic  reflection  profiles,  particularly  as  to  the 
degree  of  lateral  continuity  in  the  seabed  physical  properties. 
Logical  extension  of  the  technique  to  higher  frequencies  will 
provide  higher  resolution  for  assessing  the  more  surficial  seabed 
characteristics . 
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I NTRODUCT I ON 


When  we  attempt  to  interpret  continuous  seismic  reflection 
profiles  it  is  not  easy  to  ignore  the  desire  and  need  to  classify 
in  some  way  the  differences  in  physical  properties  of  the  sediment 
and  rock  that  cause  the  reflection  horizons  that  are  observed. 

Much  of  such  classification  is  based  upon  critical  judgment.  A 
significant  part  is  based  upon  determinations  of  the  interval  and 
critical-refraction  velocities  of  the  sediments  and  rock  derived 
from  wide-angle  reflection-refraction  profiles  obtained  with  sono- 
buoys  (LePichon  et  al.,  1968),  upon  sampled  outcrops,  and  upon  the 
cores  retrieved  and  the  down  hole  measurements  made  at  drill-hole 
sites,  particularly  those  of  the  Deep  Sea  Drilling  Program  (DSDP) . 

Sonobuoys  and  drill-holes  are  site-or  station-type  obser- 
vations so  it  becomes  a question  as  to  how  far  one  can  extrapolate 
the  classifications  at  such  sites  to  other  locations  primarily  on 
the  basis  of  the  continuity  of  reflection  horizons.  Even  if  good 
continuity  is  observed  there  are  often  large  apparent  changes  in 
reflected  energy  that  cause  one  to  wonder  how  constant  or  variable 
the  contrast  in  materials  may  be  along  a given  horizon  traced  from 
one  location  to  another. 

In  this  paper  we  describe  how  we  have  evolved  and  used  a 
relatively  straight-forward  method  for  systematically  deriving 
estimates  of  the  acoustic  impedance  of  the  seabed  materials  and 
the  impedance  variation  with  depth  and  location  using  low  fre- 
quency normal-incidence  reflections.  ("Seabed"is  used  to  imply 
everything  below  the  sea  floor.)  We  have  sought  such  estimates 
not  only  to  improve  our  ability  to  extrapolate  classifications  be- 
tween sampling  sites,  but  also  these  measurements  like  those  of 
velocity  provide  an  added  dimension  which  must  be  justified  within 
the  context  of  the  physical  and  geological  interpretations  of  the 
reflection  profiles.  It  is  readily  understood  that  knowledge  of 
velocity-depth  variations  is  essential  to  the  determination  of 
sediment  thicknesses,  to  the  quantitative  assessment  of  sediment 
volume  and  to  geologic  inferences  that  relate  velocity  to  type  of 
material,  age  and  consolidation.  Estimates  of  acoustic  impedance 
complement,  corroborate  or  deny,  the  same  suite  of  critical  judg- 
ments and  classifications  giving  better  definition  to  our  in- 
ferences. 
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To  demonstrate  the  use  of  this  method  we  have  examined 
seismic  reflection  energetics  at  six  locations  (sonobuoy  stations 
#3,  #4,  #8,  #14,  #15  and  #17)  in  or  near  the  environs  of  the 
Somali  Basin  (Fig.  1) . We  chose  to  make  this  study  of  the  Somali 
because  the  wide  extent  of  somewhat  similar  sedimentary  structure 
of  varying  thicknesses  above  oceanic  basement  permits  relatively 
simple  comparisons  of  reflection  energetics,  and  because  cores 
from  DSDP  #235,  near  station  #14,  cover  most  of  the  sedimentary 
structure  down  to  basement  (Fisher,  Bunce  et  al.  1974).  Wide- 
angle  and  normal  incidence  profile  data,  tape  recorded  at  each 
sonobuoy  station,  are  used  in  the  analyses,  wide-angle  for  velocity 
structure,  and  normal  incidence  for  reflectivity  measurements  from 
which  impedance  estimates  are  derived,  and  from  which  the  continuity 
of  the  contrasts  in  seabed  materials  is  evaluated. 


BACKGROUND 


Mapping  the  distribution,  thickness  and  structure  of  the 
sedimentary  accumulations  and  underlying  crustal  structures  below 
the  seafloor  are  important  scientific  objectives.  The  principle 
techniques  involved,  normal- incidence  and  wide-angle  reflection 
profiling  and  refraction  profiling,  are  each  based  primarily  upon 
acoustic  travel  time  measurements.  The  propagation  velocities 
determined  by  the  latter  two  techniques  are  used  with  the  first  to 
determine  structure  thicknesses  and  shapes,  and  to  characterize 
the  materials  involved  by  their  propagation  velocities. 

Of  the  three  techniques,  normal-incidence  profiling  has  been 
the  most  extensively  employed,  probably  because  it  produces  a 
readily  assayable  graphic  record.  The  shot-to-shot  continuity  and 
coherence  in  travel  time  of  the  reflection  wavelets  along  a pro- 
file has  made  possible  the  delineation  of  the  reflection  horizons 
that  are  inferred  to  be  representative  of  geological  structure 
occurring  beneath  the  ocean  floor. 

Thus,  there  has  evolved  a sizable  literature  on  the  interp- 
retation and  classification  of  the  observed  reflection  horizons, 
on  speculations  as  to  the  causes  of  the  observed  structures,  and 
the  correlation  of  the  reflection  horizons  to  particular  geological 
structures  revealed  by  coring  and  deep  sea  drilling.  The  effective- 
ness of  the  seismic  reflection  technique  was  demonstrated  (Hersey 
& Ewing,  1949;  Officer,  1955;  Shor,  1959;  Hersey,  1963;  Ewing,  J. , 
1963),  applied  to  regional  studies  (e.g.  Ewing,  J.,  et  al,  1962; 


Ewing  & Ewing,  1962;  Moore  and  Curray,  1963;  Bunce  and  Hersey, 
1966;  Emery  and  Uchupi,  1965;  Rona  and  Clay,  1967;  Knott  and 
Hoskins,  1968)  and  results  synthesized  (Ewing,  et  al,  1964; 
Talwani,  et  al,  1972).  Nafe  & Drake  (1957),  Breslau  (1966), 
Hamilton  (1970,71,76),  Houtz,  Ewing  & LePichon  (1968),  Whitmarsh 
(1972)  and  others  provide  a wide  foundation  of  insight  on  the 
velocity  and  density  structure  of  oceanic  crust,  on  seafloor  re- 
flectivity, on  the  relationships  between  reflectivity  and  the 
surficial  physical  properties  of  the  seafloor,  and  on  the  com- 
pressional  wave  velocity  characteristics  of  the  sediments.  Comp- 
lementary studies  of  the  seismic  energy  transmitted  via  reflected 
and  refracted  paths  have  also  a long  history  (e.g.  Knott,  1899; 
Nafe,  1957;  Gupta,  1965,  Hermont,  1967;  Claebout,  1968;  Borcherdt, 
1973. 


But  to  our  knowledge  the  reflection  energetics  associated 
with  the  now-conventional  continuous  reflection  profile  has  lacked 
the  systematic  examination  that  has,  for  instance, been  given  echo 
sounding  of  the  seafloor  (e.g.  Breslau,  ibid,  and  Tyce,  1976).  It 
is  convenient  now  to  use  digital  techniques  to  display  and  measure 
the  distribution  of  reflected  energy  with  travel  time.  Thus  we 
can  more  readily  take  advantage  of  these  measurements  to  increase 
our  knowledge  of  the  physical  changes  in  seabed  characteristics 
evidenced  in  the  differences  in  reflection  energetics  as  functions 
of  depth  and  region. 


The  shot-by-shot  pressure  time  curves  comprising  the  re- 
flected wavelet  sequences  from  the  seabed  are  the  primary  inputs 
to  the  analyses  (Fig.  2),  but  also  a knowledge  of  the  reflection 
structure  along  a traverse  at  the  location  is  required  so  that 
phenomena  such  as  side  echoes  and  basement  relief  can  be  identified 
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The  data  manipulation  concepts  that  we  have  found  necessary 
for  such  a study  and  the  mathematical  justification  for  the 
physical  relationships  between  the  measured  reflection  response 
and  impedance  estimates  of  the  seabed  are  described  first  in  this 
section.  The  analyses  are  applied  to  traverse  reflection  profiles, 
and  methods  for  mapping  seabed  characteristics  are  considered. 


ressure-time  curves  and  reflection  profiles 
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and  avoided  if  possible  (Fig.  3)  . 


The  Cumulative  Energy  Curve 

Our  strategy  for  displaying,  examining  and  measuring  the 
changes  in  reflectivity  with  depth  and  region  is  to  use  the  rate, 
in  respect  to  travel  time,  at  which  reflected  energy  is  returned. 
This  is  represented  by  the  cumulative  reflected  energy  curve. 

Figure  2.  Such  curves  provide  a quickly  derived  definitive  and 
accurate  assessment  of  the  relative  strength  or  magnitude  of  the 
successive  reflection  events.  The  use  of  these  curves,  generated 
by  squaring  and  integrating  the  pressure  time  curves,  was  intro- 
duced in  underwater  acoustic  studies  at  Woods  Hole  about  two 
decades  ago  by  C.  B.  Officer  and  F.  T.  Dietz  (personal  association) 
and  Baxter  (1958,  60) . Breslau  (ibid)  took  advantage  of  their  use 
for  seafloor  reflectivity  measurements,  as  have  Knott  et  al  (1972) 
and  Baggeroer  et  al  (1973)  in  preliminary  assessments  of  the  seis- 
mic reflected  energy  budget  of  the  seabed. 

A cursory  examination  of  cumulative  energy  curves  from 
various  locations  demonstrates  that  the  time  rate  of  energy  re- 
ception varies  markedly  from  location  to  location.  Fig.  4.  Of 
interest  is  the  distribution  in  time  and  in  magnitude  of  the  series 
of  steps  and  ramps  corresponding  to  reflected  energy  arrivals  from 
the  seafloor  and  deeper  interfaces.  At  some  locations  the  re- 
flection from  the  seafloor  is  the  largest  return,  at  other  locations 
a deeper  interface  may  give  the  strongest  return.  The  size  of  the 
steps  or  the  rate  of  cumulation  decreases  with  travel  time  until 
finally  no  more  energy  received  from  depth  is  discernible  above 
noise . 

Noise 


Reflection  signals  are  mixed  with  a continuum  of  background 
noise  due  primarily  to  the  self-noise  of  the  towed  receiving  hydro- 
phone array.  Integration  of  the  noise  alone  produces  a near  linear 
ramp.  To  remove  this  ramp  on  which  the  reflected  signals  would 
appear,  the  noise  energy  is  measured  per  unit  time  before  the  sea- 
floor return  (b  in  Fig.  2),  and  this  amount  is  subtracted  per  unit 
time  during  the  integration  of  seabed  returns.  Thus  the  signal-to- 
noise  ratio  during  the  seabed  returns  is  in  rea  lity  unity  or  0 dB 
whenever  the  trace  is  horizontal  although  the  ratio  may  exceed 
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FIG.  3 (a)  Reflection  profiler  recordings  at  the  six  stations. 

(b)  Profiler  recording  showing  continuity  in  sediment  reflections  between  statior.s 
#3  and  #4,  adjacent  to  Chain  Ridge  basement  outcrops. 


Cumulative  energy  traces  from  several  sites.  (Maximum  deflection 
normalized  to  a standard  unit  deflection). 
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20  dB  for  the  sea  floor  reflection.  Obviously  the  lower  the 
noise  the  better  the  determinations. 


Normalization  of  trace 

Note  that  the  cumulative  energy  curves  (integrands,  for 
brevity)  that  are  displayed  for  simple  comparison  in  Figure  4 
are  normalized  to  a common  terminal  deflection,  whereas  in  Figure 
2 the  normalization  point  of  the  integrand  is  mid-trace  and  the 
magnitude  (relative  to  the  quiescent  trace)  of  cumulative  energy 
at  that  point  is  given  (c  in  Fig.  2). 

In  general , normalization  provides  a means  to  control  the 
size  of  recordings,  and  combined  with  measurements  of  the  ampli- 
tude or  energy  at  the  normalization  point  (norm  point)  it  estab- 
lishes scale.  Most  important  to  this  study,  normalization  rela- 
tive to  a common  basis  permits  comparisons  of  the  relative  rates 
of  change  in  magnitude  of  the  cumulative  reflected  energy  from 
one  shot  to  that  of  another  in  a situation  in  which  the  actual 
energy  of  each  shot  is  unknown.  Within  a given  bandwidth  the 
relative  energetic  response  of  the  seabed  at  a given  location  can 
be  expected  to  be  reasonably  repeatable  from  shot  to  shot  even  if 
the  absolute  shot-to-shot  energy  levels  are  unknown  but  within 
practical  limits.  Thus  similar  shaped  integrands  are  expected  from 
shot  to  shot  provided  the  curves  are  similarly  normalized  (Fig.  5) . 
The  selected  position  of  the  norm  point  in  travel  time  establishes 
a unit  reference  deflection  representing  the  energy  cumulated  from 
the  onset  of  the  seafloor  returns  to  the  reception  of  some 
particular  reference  reflection  horizon,  such  as,  oceanic  basement, 
an  intermediate  sediment  horizon,  the  sea  floor,  or  between  hori- 
zons at  a quiescent  point  in  the  traces.  Like  most  acoustic  data, 
shot-to-shot  variations  are  expected  in  the  energy  cumulated  up  to 
the  norm  point.  These  variations  tend  to  become  greater  with  in- 
creasing travel  time  because  of  the  cumulative  effect  of  early 
variations  upon  those  arrivals  that  occur  later,  and  because  the 
roughness  of  reflection  horizons  commonly  increases  with  depth. 

Thus  the  integrand  that  typifies  the  response  at  a given  location 
is  derived  by  smoothing  the  data. 

Two  methods  were  examined  for  deriving  a representative 
average-shape  integrand,  1)  an  average  of  similarly  normalized 
integrands,  2)  the  integrand  of  the  average  of  several  pressure- 
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Graphic  comparison  of  the  response  of  ten  shots  (stations 
i“3a  and  #8)  from  which  an  average-shape  integrand  is  com- 
puted. The  traces  are  aligned  in  travel  time  on  the  sea- 
floor reflection  and  have  a common  deflection  at  the  norm 
point.  The  average-shape  integrand  is  the  representative 
cumulative  energy  curve  for  each  station.  Data  for  station 
#3  was  taken  from  the  beginning  of  profile  (Fig.  3). 


time  curves.  The  first  option  produces  the  better  representation 
for  the  data  at  hand  and  is  used  in  this  study.  Figure  5.  Al- 
though an  average  of  amplitude  traces  improves  the  signal-to-noise 
ratio  for  coherent  signals  such  as  the  reflections  from  the 
horizontal-lying  sediments,  the  fluctuation  in  the  travel  times 
of  signals  from  a rough  interface  such  as  basement  causes  these 
reflection  signals  to  add  incoherently  producing  less  than  actual 
amplitudes  for  these  signals.  if  the  data  that  are  averaged  were 
collected  in  one  spot-location  rather  than  along  a short  section 
of  a traverse,  the  second  option  might  be  superior. 


Frequency  effects 

There  are  frequency  dependent  effects,  and  in  order  to  re- 
duce their  significance,  all  data  are  filtered  with  a low  fre- 
quency, 30  to  60  Hz,  bandpass.  The  wavelengths  range  from  about 
20  to  90  meters  in  the  sediments  depending  on  the  propagation 
velocity.  The  reflection  of  a signal  having  at  least  an  octave 
bandwidth  and  long  wave  lengths  is  less  subject  to  major  amplitude 
changes  due  to  interference  effects  resulting  from  closely  spaced 
reflection  horizons  than  narrower  band  and  shorter  wavelength  sig- 
nals. Absorption  is  low; and  incoherent  backscattering  from  sedi- 
ment horizons  is  less  than  at  higher  frequencies  because  the 
magnitude  of  the  horizon  roughness  is  a smaller  fraction  of  the 
longer  wavelengths. 


Integrand  comparisons 

In  order  to  compare  the  integrands  other  than  in  the 
cursory  manner  of  Figure  4,  let  us  divide  the  integrand  de- 
flection into  equal  increments  of  cumulative  energy  and  determine 
in  descriptive  terms  how  the  increments  may  be  related  to  the 
characteristics  of  the  seabed.  in  the  discussion  of  normalization 
we  have  pointed  out  the  advantage  of  normalization  relative  to  some 
reference  horizon.  If  for  instance  the  reference  horizon  were 
oceanic  basement,  we  can  divide  the  integrand  into  equal  deflection 
increments  relative  to  the  integrand  deflection  immediately  after 
the  reception  of  the  basement  return  and  assign  percentages  to  the 
increments  relative  to  the  "total”  energy  cumulated  to  this  point. 
Thus  we  assume  a synthetic  series  of  equal  energy  steps  and  then 
examine  ways  to  relate  the  reflected  energy  step-by-step  to  what 
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the  physical  characteristics  of  the  seabed  would  have  to  be  to 
cause  each  step. 

The  concept  of  generalizing  the  reflection  integrand  by 
representing  it  as  a series  of  steps  is  perhaps  unique.  We 
rationalize  this  concept  on  the  basis  that  the  percentage  values 
can  be  made  small  to  fit  an  integrand  having  a number  of  discrete 
arrivals,  and  this  approach  simplifies  the  mathematical  analyses. 

We  know  rather  well  that  the  acoustic  impedance  of  sea  water 
at  depth  is  about  1,6  x 10°  rayls  and  of  basaltic  oceanic  basement, 
about  16  to  18  x 10b  rayls.  Therefore  in  our  mathematical  and 
physical  analyses  we  seek  to  relate  these  equal  energy  percentile 
steps  to  changes  in  the  acoustic  impedance  between  these  boundaries. 

The  equal  steps  of  cumulative  energy  can  also  be  calibrated 
in  relation  to  the  reflectivity  of  the  sea  floor  once  reflectivity 
is  measured  or  assumed.  The  deflection  of  the  first  event  of  the 
integrands  can  then  be  made  to  represent  the  sound  power  reflection 
coefficient  assuming  the  response  of  the  seafloor  is  the  same  from 
shot-to-shot.  Essentially  this  assumes  that  energy  incident  at  the 
seafloor  is  unity  and  the  first  deflection  some  fraction  (the  re- 
flection coefficient)  of  the  incident  energy.  The  cumulative  energy 
at  each  integer  multiple  of  the  first  deflection  can  then  be 
measured  relative  to  the  coefficient  thus  again  configuring  an  in- 
tegrand of  equal  energy  steps. 

By  assuming  incident  energy  is  unity,  we  are  not  dependent 
upon  a knowledge  of  the  shot-to-shot  variation  in  seismic  source 
level,  but  do  depend  in  all  these  measurements  upon  a reasonable 
linearity  in  the  response  of  the  seabed  to  variations  in  incident 
energy.  Baggeroer  and  Hoskins  (1975)  measured  in  a 10-400  Hz  band- 
pass fluctuations  of  2 to  6 dB  in  the  source  level  of  a 300-cubic 
inch  airgun  using  a reference  hydrophone  46  meters  directly  below 
the  source.  Knott  et  al  (ibid)  observed  in  a 20-120  Hz  bandpass 
variations  of  3 dB  in  the  downward  radiated  energy  from  the  100 
kj  underwater  spark,  measured  4.5  meters  below  the  source  in  sea 
state  1 to  2 . Figure  5 shows  at  two  locations  that  the  relative 
seabed  response  is  quite  repeatable  for  ten  successive  shots. 

This  suggests  linear  response  at  least  over  the  2 dB  range  of  the 
energies  measured  at  the  norm  point  of  the  original  integrands 
from  which  this  figure  was  produced. 


Mathematical  description  of  the  reflection  model 


A quantitative  analysis  of  the  sequence  of  acoustic  re- 
flections from  the  seabed  needs  to  begin  with  a logical  and 
causal  mathematical  model.  Although  the  model  described  here 
does  not  incorporate  all  the  phenomena  involved  in  the  reflection 
and  transmission  through  a sequence  of  seabed  layers,  it  does 
provide  an  adequate  identification  of  and  an  explanation  for  most 
of  the  received  signal.  We  will  first  examine  a steady  state 
model  and  then  show  that  geophysical  realities  allow  major  simpli- 
fication of  this  model  so  that  it  is  suitable  for  the  case  of 
transient  signals.  Since  this  mathematical  model  is  based  on  a 
simplification  of  the  actual  situation,  those  geophysical  features 
that  appear  to  be  of  secondary  importance  and  that  can  be  neglected 
need  to  be  stated. 

The  acoustic  pulse  travels  from  a source  near  the  surface  to 
the  seafloor  and  returns  to  a receiver  also  near  the  surface  and 
the  source.  The  events  of  interest  here  are  the  reflections  from 
the  water-seabed  interface  and  the  reflections  from  within  the  sea- 
bed. The  transmission  path  from  the  surface  to  the  seafloor  and 
the  return  path  from  the  seafloor  is  virtually  the  same  for  all  the 
reflections  of  one  data  set;  the  only  difference  in  paths  for  these 
different  reflections  is  within  the  seabed  itself. 

Within  the!  seabed  the  sound  is  assumed  to  propagate  as  a 
plane  wave.  For  deep  sea  seismic  profiling  the  depth  of  the  water 
is  greater  than  the  depth  of  penetration  into  the  seabed.  Thus, 
the  transmission  loss  due  to  inverse-square  spreading  within  the 
seabed  is  small  and  will  be  neglected  for  the  purpose  of  this  re- 
port. 


The  incident  and  reflected  wave  paths  are  assumed  to  be  at 
normal  incidence  to  the  interfaces.  The  angles  that  do  occur  are 
so  small  that  any  corrections  due  to  Snell 1 s Law  are  much  less  than 
the  resolution  in  travel  time.  Compressional  wave  propagation  is 
assumed . 

We  assume  the  losses  from  scattering,  absorption,  and 
spreading  within  the  water  column  are  the  same  for  all  subject 
transmission  paths.  Variations  in  these  common  losses  drop  out 
by  relating  the  reflection  signals  to  one  another  in  a relative 


-14- 


manner  or  to  the  sound  incident  at  the  water-seafloor  interface 
by  means  of  seafloor  reflectivity  measurements. 

Within  the  seabed  any  losses  due  to  absorption  and  scat- 
tering are  assumed  to  be  small  in  comparison  to  the  effects  due  to 
the  reflections.  Both  absorption  and  scattering  are  frequency 
dependent  and  are  small  for  low  frequencies,  i.e.  a 30  to  60  Hz 
bandpass . 

Within  a given  layer  the  acoustic  properties  are  assumed  to 
be  homogeneous  and  isotropic,  with  discrete,  non-reactive  inter- 
faces. Further,  we  assume  that  the  seabed  responds  linearly  over 
the  range  of  shot-to-shot  variations  in  acoustic  pressure  generated 
by  the  seismic  sources. 

Implicit  in  the  interpretations  of  measurements  of  reflected 
energy  versus  travel  time  is  the  ambiguity  as  to  whether  the  re- 
flections observed  are  due  to  an  increase  or  decrease  of  impedance 
between  succeeding  layers.  An  overall  increase  in  acoustic  im- 
pedance with  depth  is  to  be  expected,  but  there  may  be  some  im- 
pedance decreases  or  "inversions"  interspersed.  At  the  present 
we  assume  that  there  are  no  inversions.  Analysis  based  on  this 
assumption  gives  an  upper  bound  but  not  a lower  bound  to  the  value  of 
acoustic  impedance  versus  depth  (Fig.  6c) . 

Several  factors  enter  into  what  produces  changes  in  acoustic 
impedance  and  thereby  a reflection.  The  most  easily  visualized 
geological  structure  is  one  made  up  of  planar  boundaries  between 
materials  of  differing  density  and/or  compressional  wave  velocity. 

In  nature,  boundaries  may  be  discrete  or  transitional,  or  more 
specifically,  discrete  or  transitional  in  respect  to  the  band- 
width and  wavelength  of  the  signal  insonifying  it  (Wolf, 1937). 

One  can  also  anticipate  that  the  reflection  wave  train  sequence 
observed  will  differ  depending  on  the  characteristics  of  the  in- 
cident waveform.  The  smoothness,  slope,  and  extensiveness  of  the 
interface  also  affect  the  manner  in  which  it  reflects  sound. 

For  convenience  in  the  development  of  the  model  we  consider 
three  separate  situations  for  steady  state  waves:  the  first,  a 

single  interface,  the  water-seafloor  reflector;  the  second,  a 
basaltic  basement  with  only  one  layer  above  it;  and  the  third  or 
general  case,  several  layers  between  water  and  basement.  The 
first  situation  is  covered  in  acoustic  texts  (e.g.  Kinsler  & Frey, 
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Application  of  the  five  20-percentile  synthetic  layer  model 

to  data. 

(a)  Cumulative  energy  curve  re  travel  times  to  20$,  ^0$, 
60$,  80$  and  100$  of  the  energy  cumulated  up  through 
reception  of  the  basement  return. 

(b)  Five  synthetic  layers  with  interfaces  corresponding 
to  each  equal  percentile  step,  labelled  with  maxi- 
mum impedance  predicted  from  model.  (*  18.0  repre- 
sents a high  average,  near  maximum  impedance  for 
basalt . ) 

(c)  Plot  of  estimated  impedance  vs  travel  time.  The 
impedance  is  predicted  to  be  between  the  solid  and 
dashed  line.  Maximum  values  (dashed-line)  are  pre- 
sented in  all  other  figures. 
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1962).  The  second  case  is  usually  treated  in  texts  (ibid  p.136) 
as  a transmission  problem.  We  will  devote  considerable  dis- 
cussion to  this  case  because  it  will  allow  us  to  make  some  cal- 
culations which  produce  some  simplifications  applicable  to  the 
third  case.  Then  the  third  case  can  be  related  to  reflection  of 
acoustic  pulses. 

For  a plane  progressive  wave  at  normal  incidence  on  a bound- 
ary between  the  water  and  the  seafloor,  a reflected  wave  and  a 
transmitted  wave  are  generated.  The  ratio  of  the  amplitude  of  the 
primary  reflected  wave,  &i  , to  the  incident  wave,  is 


<°A  + 


2,  - / 

* 


where  PpD  and  /°P>  are  the  characteristic  impedances  of  the  water 
and  the  seafloor  respectively  and  / /0ac-v  . The  sound 

power  reflection  coefficient  d,  is  the  ratio  of  the  intensity  of 
the  reflected  wave  to  the  intensity  of  the  incident  wave  or 


PpPPP-j 


Eq.  1 


For  the  layered  media,  a careful  distinction  must  be  main- 
tained between  the  reflection  process  at  an  individual  interface 
and  the  reflection  process  involving  all  the  layers.  For  an  indi- 
vidual interface  we  use  the  terms  reflectivity,  r,  and  sound  power 
reflection  coefficient,  <A  (see  Eq.  1) . The  effective  sound  power 
reflection  coefficient  for  an  entire  sequence  will  be  designated 
as  the  reflectance,  , and  it  is  defined  as  the  ratio  of  the  net 
intensity  of  the  reflected  wave  to  the  intensity  of  the  incident 
wave.  For  the  case  of  only  one  interface,  the  reflectance  is  the 
sound  power  reflection  coefficient  or^,=  ct , . For  more  than  one 
interface,  the  reflectance  depends  on  the  sound  power  reflection 
coefficients  of  all  the  interfaces. 


Next  let  us  consider  the  second  case  when  the  incident  plane 
wave  from  the  water  is  transmitted  through  a finite  layer  into  the 
third  medium.  When  the  incident  wave  arrives  at  the  first  inter- 
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face  the  wave  is  partly  reflected  and  partly  transmitted  into 
the  second  medium.  This  transmitted  wave  continues  until  it 
arrives  at  the  second  interface.  Here  it  is  partly  transmitted. 

The  reflected  wave  returns  to  the  first  interface  where  again 
partial  reflection  and  partial  transmission  into  the  water  occurs. 
This  transmitted  wave  combines  with  the  wave  initially  reflected 
in  the  water,  while  the  reflected  wave  combines  with  the  wave 
initially  transmitted,  creating  a steady  state  condition. 

The  wave  amplitudes  are  summed  with  their  proper  phase  re- 
lationship according  to  the  Principle  of  Superposition  for  linear 
waves.  The  energy  or  the  power  is  then  computed  from  the  resultant 
amplitude.  The  phase  change,  S • between  successive  reflections 
is  the  same,  2 IT  (path  difference/wavelength),  or  2k^h^where 
v = 2U/\i,  and  is  the  thickness  of  the  layer.  The  reflectance 
for  the  single  layered  seabed  in  terms  of  6 and  the  sound  power 
reflection  coefficients0^, and  o(L  is 

oL,  2 - cos  Jr 

, r Eq.  2 

/ -A  d ( - 2 <^2  Cus  S 

(Details  for  the  derivation  of  this  equation  are  given  in  Appendix 
I)  • 


The  reflectance  given  in  Eq.  2 is  very  dependent  on  the  phase 
change  S . For  a given  thickness,  the  reflectance  depends  on 
the  frequency.  In  optical  applications  and  in  most  acoustical 
situations  the  calculations  of  the  reflectance  are  made  for  a 
single  frequency  or  a narrow  band  of  frequencies.  However,  an  im- 
portant distinction  must  be  made  in  the  case  we  are  dealing  with, 
the  signal  is  broad  band  covering  at  least  one  octave.  Thus  the 
effective  reflectance  that  will  be  measured  is  the  average  over 
the  frequency  band.  For  typical  values  of  (A,  and  <X2  , the  product 
<p( t c*  2.  in  the  denominator  of  Eq.  2 is  very  much  less  than  1 

so  it  may  be  neglected.  To  simplify  Eq.  2 for  averaging  over  the 
frequency  band,  the  denominator  is  expanded  by  the  binomial  theorem 
keeping  terms  on  the  order  of  <dz  •. 

(ff  ^ £*,-/•  O' 2.  ~Z  Vdf&z  CdiJ  + $ A,  Ces  </ 

But  the  average  of  cos  S over  the  frequency  band  is  zero  while  the 
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average  of  cos  d is  1/2.  Thus  the  effective  value  of  the  re- 
flectance to  the  level  of  for  the  broad  band  signal  is 


V ot,  + ^2  ^ 


Eq.  3 


In  Appendix  II  a comparison  between  the  average  values  of  (R.2  com- 
puted by  Eq.  2 and  Eq.  3 is  made  for  values  of  ^7  / p,  C0  between 
1 and  10  (where  P*  / p„C0  = 10).  This  is  equivalent  to  a 
variety  of  layers  between  water  and  the  basaltic  basement  rock. 

The  difference  between  the  calculations  from  Eqs.  2 and  3 is  less 
than  the  expected  errors  in  the  measured  values  of  c and  the  re- 
flectance. 


1 


Equation  3 can  be  simply  derived  another  way.  The  averaging 
over  the  frequencies  implies  that  the  phase  variations  for  broad 
band  signals  is  not  critical  and  hence  we  can  deal  with  intensi- 
ties directly.  Also  by  keeping  terms  to  the  order  of  we  need 

to  consider  the  initial  reflections  and  the  first  reflections  of 
the  down  going  rays  only,  i.e.,  returned  energy  from  "multiples" 
can  be  ignored.  Thus  the  intensity  of  the  reflected  signal  will 
be  the  sum  of  the  intensity  of  the  signal  reflected  from  the  first 
surface  plus  the  signal  transmitted  through  the  first  surface  re- 
flected from  the  second  and  thence  into  the  water  again,  or 


, 1 o -f  Z"  1 o 
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where 


A 


is  the  transmission  coefficient. 


But  - 1 - U,  so 

(R±  & -h  (P2  Eq.  3(a) 


This  formulation  will  greatly  simplify  the  calculations  for  the 
general  case. 


■v  A., 


WiWNM wm 


-19- 

The  general  solution  for  the  multi-layered  seabed  is  more 
complicated.  If  the  method  outlined  for  Eq.  2 is  followed  (see 
Appendix  I),  each  reflection  generates  an  infinite  series  to  be 
summed,  quickly  producing  an  unmanageable  form.  The  other 
approach  (see  Fowles,  1975)  is  to  set  up  a transfer  matrix  for 
each  layer.  Then  the  product  of  the  individual  matrices  yields  an 
overall  transfer  matrix  from  which  the  reflectance  can  be  computed — 
again  a difficult  task  for  the  general  case.  However,  if  the 
simplifications  that  allowed  us  to  calculate  Eq.  3(a)  are  used, 
the  general  equation  for  the  reflectance  is  readily  computed.  The 
details  are  given  in  Appendix  III.  For  four  layers  with  five 
interfaces  the  total  reflectance  can  be  written  as 


Eq.  4 


where  is  the  reflection  coefficient  between  successive  layers 
(Eq.  1)  and  T is  the  corresponding  transmission  coefficient.  The 
subscripts  designate  the  successive  interfaces. 

All  of  the  preceding  analysis  is  based  on  a steady  state 
solution  while  the  experimental  signals  are  broad  band  pulses. 
However,  the  conditions  which  permitted  us  to  write  Eq.  4 are 
equally  applicable  for  broad  band  pulses.  Averaging  over  a fre- 
quency band  essentially  decouples  the  phase  variations  for  the 
steady  state  calculation;  the  same  should  be  true  for  pulse  oper- 
ation. By  keeping  terms  to  the  order  of  only,  the  multiple  re- 
flections are  assumed  to  produce  a negligible  contribution  to  the 
received  signal  and  therefore  the  reflection  from  each  interface 
is  independent  of  the  other  reflections.  The  analysis  used  for 
Eq.  4 should  apply  equally  well  for  pulses.  In  other  words  the 
total  reflected  signal  can  be  considered  to  be  the  sum  of  a sequence 
of  independent  reflections  from  successive  interfaces  (Appendix 
III);  hence  the  use  of  the  cumulative  energy  curve. 

For  computational  purposes  we  have  initially  structured  our 
model  to  represent  the  acoustic  impedance  contrasts  that  would 
cause  five  equal  (20%)  energy  increments  (Fig.  6).  The  intervals 
then  in  effect  become  synthetic  equivalent  acoustic  impedance 
steps,  or  generalized  "horizons"  between  generalized  "layers". 

The  increments  in  travel  time  and  hence  the  apparent  thickness  of 
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the  generalized  layers  depend  on  the  shape  of  the  observed 
cumulative  energy  curve,  Figure  6a.  The  impedance  modelled  in 
our  study  area  is  constrained  by  relatively  well  known  bound- 
aries, at  the  top  of  the  sediment  column  by  the  impedance  of 
seawater  and  at  the  bottom  by  that  of  basaltic  basement.  Energy 
percentiles  on  the  cumulative  energy  curve  range  therefore  from 
zero  at  the  onset  of  sediment  returns  to  100%  after  reception  of 
the  basement  return. 


The  model  is  applied  in  two  ways,  and  four  conditions  con- 
trol its  effectiveness  for  estimating  the  acoustic  impedance 
structure  with  depth.  In  the  first,  the  model  predicts  a series 
of  impedance  value  steps  for  each  20%  step  in  cumulative  energy. 

In  the  second  it  predicts  a series  of  impedance  steps  for  integer 
multiples  of  the  sound  power  reflection  coefficient,  , in 

terms  of  reflectance.  Effectiveness  is  controlled  by: 

1.  Where  there  are  both  a return  from  basaltic  base- 
ment and  a determination  of  seafloor  reflectivity, 
the  model  can  be  applied  both  ways  for  comparisons. 

2.  Where  there  is  a basement  return,  but  no  reflect- 
ivity measurement,  the  first  method  is  applied. 

3.  Where  there  is  no  discernible  basement  return,  but 
reflectivity  is  measured,  the  second  method  is  applied. 
The  first  may  be  used  as  an  approximation,  but  a lesser 
impedance  value  than  that  of  basalt  must  be  assumed 
for  the  terminal  value  of  the  cumulative  energy  trace. 

4.  Where  there  are  neither  a basement  return  nor  a re- 
flectivity measurement,  relative  impedance  variations 
can  only  be  estimated  on  the  basis  of  geological  or 
reflectivity  and  reflectance  assumptions. 

The  reflection  model  can  be  used  to  estimate  the  impedance 
structure  and  to  establish  a base  from  which  to  measure  anomalous 
changes  witnessed  in  the  reflection  energetics.  Even  if  the  model 
is  not  precise  it  can  thus  provide  a basis  for  comparisons  between 
the  energetics  observed  in  different  areas. 

Impedance  re  reflectance 


In  order  to  work  with  equal  steps  of  cumulative  energy  the 


reflected  contributions  from  each  interface  defined  in  Eq.  4 
must  be  made  equal (Figure  III-l,  Appendix  III).  in  Eq.  4 the 
second  term  ( T,  ) is  made  equal  to  the  first  ( ) , the  third 

( rtf  ) equal  to  the  second  ( ) and  so  on.  In  this  way  by 

starting  with  an  appropriate  £>(,  or  sound  power  reflection  co- 
efficient for  the  seafloor  we  can  solve  for  the  successive  's. 

Then  the  impedance  contrasts  and  subsequently  the  impedances  can 
be  derived  from  Eq.  1.  A family  of  curves  is  then  generated  by 
substituting  a range  of  cK  ' s which  gives  the  acoustic  impedance 
of  each  synthetic  equivalent  layer  and  the  number  of  synthetic 
layers  to  be  expected  between  the  seawater  and  basaltic  basement 
(Fig.  7 and  6) . 

The  logic  is  reversed  to  relate  impedance  to  the  percentage 
of  total  reflected  energy.  The  impedance  of  the  5th  20-percentile 
synthetic  equivalent  layer  has  then  to  be  considered  in  terms  of 
the  maximum  it  will  reach  in  that  interval,  terminating  at  the 
bottom  of  the  layer  with  an  impedance  of  basalt  (dashed  line  in 
Fig.  6c) . Further,  the  reflectance  of  the  first  layer  is  not 
equivalent  to  the  sound  power  reflection  coefficient  of  the  sea- 
floor as  above,  but  rather  a synthetic  cK  for  the  entire  first 
synthetic  equivalent  layer.  In  this  case  the  impedance  of  each 
synthetic  layer  is  read  off  that  ordinate  in  Figure  7 where  the 
line  signifying  the  5th  layer  crosses  an  appropriate  terminal  value 
of  impedance  for  basalt  (16  to  18  x 10b  rayls) . 

Impedance  structure  can  therefore  be  estimated  rather  easily 
from  much  of  the  conventional  deep-sea  seismic  profile  data.  Base- 
ment returns  are  commonly  received,  and  the  extent  of  the  oceans 
where  water-column  multiples  do  not  interfere  is  great. 

Use  of  the  seafloor  return  as  a reference  horizon  has  the 
advantage  of  an  omni-present  reference,  but  the  added  compli- 
cations associated  with  seafloor  reflectivity  measurements  are 
necessarily  involved.  Reflectance  based  upon  seafloor  reflecti- 
vity is  particularly  useful  lor  assessing  the  reflected  energy 
budget,  because  the  deflection  at  any  point  of  the  integrand  scaled 
to  the  sound  power  reflection  coefficient  equals  the  reflectance 
at  that  depth  or  travel  time.  With  such  calibrated  curves  we  can 
refer  to  the  "sediment  reflectance"  (reflectance  to  the  base  of 
the  sediments)  or  the  "basement  reflectance"  (total  reflectance  to  just 
after  reception  of  the  basement  return) . 
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lace  causing  reflections  of  equal  energy  as  observed  in  the  water  layer.  Plotted 
as  a function  of  the  sound  power  reflection  coefficient  and  reflectivity  co- 
efficient of  the  "first"  seabed  interface.  When  the  20-percentile  method  is  used 
the  "first  layer"  extends  from  the  seafloor  down  to  t20-  The  impedance  of  this 
first  layer  may  or  may  not  be  equal  to  that  of  the  seafloor. 
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Measurement  uncertainties 

In  each  approach,  either  via  reflectance  or  percentage,  it 
is  necessary  to  consider  and  define  what  constitutes  the  return 
from  the  reference  interface.  This  need  is  best  demonstrated  by 
the  reflectance  measurements  based  on  seafloor  reflectivity.  if 
the  sound  power  reflection  coefficient  is  already  known,  the 
question  is,  at  exactly  what  point  on  the  integrand  does  the  de- 
flection equal  the  coefficient?  We  define  this  point  in  time 
from  the  onset  of  returns  to  be  equal  to  the  "effective  pulse 
length"  of  the  seismic  source.  If  this  time  is  picked  too  early 
or  too  late,  the  deflection  will  be  too  low  or  too  high  and  this 
will  cause  subsequent  scaling  errors  in  reflectance.  Reasonable 
comparisons  of  reflectance  can  be  made  if  such  an  error  is  con- 
sistent because  the  error  approximates  a coefficient  applied  to 
each  reflectance  that  is  compared.  Estimates  of  impedance,  how- 
ever, will  be  greatly  distorted  (Fig.  8).  Concerning  basement 
as  a reference,  it  is  unusual  when  the  relief  in  travel  time  along 
a profile  is  not  several  times  the  period  of  a basement  return,  but 
at  a station  a sequence  of  shots  can  usually  be  selected  in  which 
the  migration  of  the  basement  return  is  minimal. 


The  "effective  pulse  length"  includes  the  effects  of  the 
multi-path  sea-surface  reflections  at  the  source  and  receiver 
upon  the  characteristic  pulse  length  of  the  source.  We  have 
estimated  the  effective  pulse  length  for  these  data  to  be  50  to 
80  milliseconds,  and  when  dealing  with  the  basement  return  we  have 
extended  the  period  to  100  ms  to  allow  for  some  travel  time  vari- 
ation. The  source  in  this  case  was  the  combined  output  of  a 100 
kJ  underwater  spark  and  a 10-cubic  inch  airgun  synchronized  to  the 
first  positive  pressure  pulse  while  each  was  towed  at  about  the 
same  depth  (8  to  9 meters) . A 50  millisecond  effective  pulse 
length  is  appropriate  for  the  spark  and  the  towing  geometry  of  the 
source  and  receiver.  The  effective  pulse  length  of  the  airgun  is 
not  so  easily  defined  because  of  the  many  oscillations  of  the 
bubble  pulse  train  (period  of  one  oscillation  is  40  ms  for  this 
size  gun).  In  this  bubble  pulse  train  the  acoustic  pressure  of 
the  first  bubble  pulse  is  roughly  6 dB  less  than  the  initial  pulse, 
and  at  this  towing  depth  destructive  interference  caused  by  the 
surface  reflection  further  reduces  the  contributions  from  all  the 
airgun  bubble  pulses.  Hence  we  have  estimated  that  the  oscillations 
are  reduced  to  a relatively  insignificant  level  by  the  end  of  a 
80  millisecond  period. 


ESTIMATED  IMPEDANCE 
ray/s  (kg/m2ssc) 


of  the  model.  A pulse  length  of  80  ms  and  a sound  power  reflection  co- 
efficient (a)  of  0.032  from  the  30  to  60  Hz  reflectivity  measurements 
combine  to  produce  this  "best"  fit. 

The  effect  of  pulse  length  on  the  impedance  estimate. 

The  effect  of  variation  in  the  determination  of  a on  the  impedance  estimate. 
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Seafloor  reflectivity  and  signal-to-noise  ratio 

The  most  straightforward  way  to  determine  the  reflectivity 
or  the  reflectance  is  to  compare  the  initial  or  total  returns 
with  the  incident  signal,  but  this  requires  measurement  of  the 
incident  waveform.  Since  this  was  not  done  we  are  obliged  to 
determine  reflectivity  indirectly.  A credible  measurement  of  the 
reflectivity  coefficient  is  obtained  by  determining  the  ratio  of 
the  amplitude  of  the  initial  event  in  the  first  water-column 
multiple  of  the  seafloor  reflection  to  the  amplitude  of  the  initial 
event  in  the  directly  received  seafloor  reflection,  allowing  for 
a factor  of  two  for  transmission  loss  and  assuming  sea  surface  re- 
flectivity to  be  100%. 

With  a favorable  signal-to-noise  ratio,  reflectivity  could 
be  so  determined  by  computer  from  the  returns  of  each  shot,  but 
in  the  suoject  data  the  signal-to-noise  ratio  of  the  initial  ar- 
rival of  the  first  water  column  multiple  is  typically  1 dB.  If 
the  signal-to-noise  is  less,  the  deflection  of  the  amplitude  trace 
can  only  be  estimated  and  the  magnitude  of  the  integrand  becomes 
indeterminable  (Appendix  V) . To  improve  reflectivity  measurements 
particularly  in  the  30-60  Hz  bandpass  we  increase  the  signal-to- 
noise  ratio  4 to  6 dB  by  averaging  an  appropriate  number  of  ampli- 
tude traces  and  processing  the  resulting  waveforms. 

Since  the  frequency  spectrums  of  the  signal  and  the  noise  are 
different,  the  signal-to-noise  ratio  is  better  in  some  frequency 
bands  than  others.  For  these  data  the  240-480  Hz  bandwidth  appears 
optimum.  Reflectivity  measurements  were  determined  in  this  band, 
in  the  30  to  60  Hz  band  and  other  octave  bands. 


Reflectance  along  reflection  profile  traverses 

The  continuity  or  transitions  in  physical  character  of  the 
reflection  horizons  evidenced  by  changes  in  reflectance  along 
tracks  crossing  different  geological  provinces  are  of  particular 
interest.  To  portray  these  changes  we  construct  an  overlay  to 
the  familiar  reflection  profile  record,  consisting  of  contours  of 
the  reflection  integrand  of  each  successive  shot  (Figs.  9 and  10) . 
The  integrands,  and  thus  the  contours,  are  normalized  as  before  to 
a reference  horizon  or  reference  travel  time  and  presented  in  terms 
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of  percentage,  or  reflectance  if  the  seafloor  reflectivity  is 
known  or  estimated.  Instead  of  determining  an  average-shape 
integrand  as  for  a station,  the  deflection  values  of  the  suc- 
cessive integrands  are  smoothed  along  the  traverse  using  a run- 
ning average  algorithm.  Although  one  can  easily  infer  estimated 
impedance  from  these  profiles,  a future  step  will  be  to  present 
contours  of  estimated  impedance. 

Both  the  percentage  and  reflectance  contours  of  data  that 
have  not  been  smoothed  (Fig.  10  b and  d)  show  significant  vari- 
ations in  travel  time.  Smoothing  the  data  serves  to  identify 
the  relative  proportions  of  short  and  long  term  fluctuations.  The 
fluctuations  in  travel  time  of  the  contours  in  the  cumulative  energy 
curve  in  Figure  10  (b)  appear  to  be  due  primarily  to  variability  in 
the  strength  of  the  return  from  the  reference  horizon  at  inter- 
mediate depth.  in  Figure  10(d)  the  fluctuations  in  travel  time 
are  apparently  due  to  localized  variability  in  seafloor  reflecti- 
vity. 


Possible  criteria  for  mapping  seabed  characteristics 

We  propose  a set  of  measurements  in  addition  to  reflectivity 
and  reflectance  which  may  be  useful  for  mapping  the  characteristics 
of  the  integrands  from  different  physiographic  regions.  These  are 
objective  rather  than  interpretative  observations  that  provide  a 
relatively  quick  and  systematic  way  for  making  areal  comparisons. 


Seabed  time  constants:  The  cumulative  energy  curves  have  a 

shape  that  is  somewhat  similar  to  a charging  curve  for  a capacitor. 
One  characterization  associated  with  such  curves  is  the  "time  con- 
stant", the  time  to  rise  from  the  onset  of  the  first  seafloor  arri- 
val to  some  specified  fraction  of  its  terminal  value.  For  a 
generalized  time  constant  characterization  we  suggest  two  time  in- 
tervals, t^Q  (time  period  for  20%  rise)  as  a measure  of  the  sea- 
floor response,  and  tg0  (80%  rise)  as  a measure  of  response  time 
of  the  entire  return.  These  time  constants  are  frequency  depen- 
dent, hence  the  passband  must  be  specified.  The  time  constants 
at  low  frequencies  are  usually  longer  than  at  higher  frequency 
passbands . 

For  reflection  sequences  including  a strong  but  variable  re- 


turn  such  as  from  basaltic  basement,  an  average  integrand  should 
be  used.  The  contoured  profiles  are  helpful  in  choosing  a 
representative  group  of  reflection  sequences  for  averaging  The 
time  constant  is  also  affected  by  the  relative  proportion  of  the 
total  return  from  the  sedimentary  and  from  the  crystalline  rocks 
(Appendix  IV) . Should  one  wish  to  consider  just  the  sedimentary 
sequence,  the  time  constants  (as  a function  of  frequency)  of  the 
sedimentary  portion  of  the  integrand  may  be  more  meaningful  than 
the  entire  return. 

Shape  of  the  energy  integrands:  The  general  shape  of  the 

cumulative  energy  curve  is  another  feature  for  characterization. 
Two  extremes  of  this  curve  are  a "tread"  function  and  a "step" 
function.  In  the  first  case  very  little  energy  is  received  fol- 
lowed by  a terminal  brust  of  energy.  In  the  other  extreme  the 
energy  rises  initially  very  fast  to  almost  80%  and  then  gradually 
approaches  80%  at  t^Q  giving  a single  step-like  trace.  An  inter- 
mediate case  is  a simple  ramp  function,  the  cumulative  energy  in- 
creases linearly  with  time  giving  a straight  line  from  zero  energy 

at  t = 0 at  the  onset  of  arrivals  to  80%  at  t . 

oO 

To  obtain  a measure  of  this  shape,  we  define  the  term 
"AreagQ"  which  is  the  ratio  of  the  area  under  the  energy-time 
curve  divided  by  the  area  under  the  simple  ramp  function  to  the 
time  tRn.  For  computational  purposes  we  use  the  trapezoidal  ap- 
proximation obtained  by  dividing  the  time  into  intervals  cor- 
responding to  the  energy  levels  of  20,  40,  60  and  80  percent. 

This  application  yields  the  equation 

Area80  = 1-75  ~ ^t20+  fc40  + t60^'/(2  fc80^ 

With  this  definition  the  values  for  Area80  should  range 
from  0 for  a tread  function,  passing  through  1 for  the  ramp 
function,  to  2 for  the  step  function.  This  approximation  actu- 
ally yields  values  between  .25  and  1.75.  This  discrepancy  is 
not  serious  since  the  values  of  AreagQ  are  for  classification  and 
comparison  purposes.  So  long  as  the  method  is  consistently  used, 
the  gross  comparisons  can  be  useful.  Figure  11  displays  the  re- 
lationship between  estimated  acoustic  impedance  and  Area^. 


Data  handling 


The  normal  incidence  reflection  data  are  taken  from  analog 
magnetic  tape  recordings  collected  in  the  field  at  sonobuoy 
stations.  The  longest  taped  profile  is  therefore  about  20  km. 

The  shipboard  receiving  and  tape  recording  systems  that  were  used 
are  linear  over  a dynamic  range  of  about  40  dB  over  a bandwidth  of 
10  Hz  to  3 kHz. 

Data  from  the  tapes  are  digitized  using  two  methods  to  con- 
trol the  digitizing  rate.  In  the  first  we  expect  from  timing 
calibrations  that  the  time  fluctuations  of  the  tape  drive  will  be 
minimally  cumulative,  and  will  not  exceed  the  order  of  0.01  second 
in  the  period  of  the  6 to  7 second  reflection  travel  time.  In  this 
case  a separate  clock  drives  the  digitizer  at  2000  Hz.  This  pro- 
cedure has  appeared  adequate  for  examinations  of  the  reflected 
energy  curves  from  the  first  seafloor  reflections,  but  becomes 
tenuous  when  exact  timing  of  the  first  water  column  multiple  of 
the  seafloor  return  is  sought  for  reflectivity  measurements.  in 
the  second  method  the  digitizing  rate  is  controlled  by  a signal 
(960  Hz)  derived  from  the  original  time  base  of  the  experiment (60Hz ) 
that  was  tape  recorded  with  the  data. 

Data  digitized  by  either  method  are  analyzed  via  our  programs 
using  a 32,000  word  Hewlett-Packard  2100  computer  with  two  1.2 
million-word  discs,  two  9-track  magnetic  tape  drives,  a CRT  termi- 
nal, an  incremental  plotter  and  line  printer  (Prada,  1976).  To 
assure  the  validity  of  the  energy  measurements,  the  program  out- 
puts were  calibrated  by  using  digitized  known  signal  inputs.  Band- 
pass filtering  was  controlled  by  calibrated  analog  (Allison,  Model 
2-AR)  filters  because  the  symmetry  of  their  response  was  better  than 
we  were  able  to  produce  by  simple  time-domain  digital  filters. 


RESULTS  AND  DISCUSSION 


The  Data 


We  have  examined  tape  recorded  normal-incidence  reflection 
data  taken  concurrently  with  wide  angle  sonobuoy  profiles  from  six 
sites  (#3,  #4,  #8,  #14,  #15,  #17,  Fig.  1).  Two  adjacent  stations 
#3  and  #4  provide  a means  to  check  the  repeatability  of  results  and 
one  station,  #14,  is  only  30  km  from  DSDP  site  #235.  On  the  original 


reflection  profile  record  there  is  good  evidence  for  the  con- 
tinuity of  the  reflection  horizons  between  the  drill  hole  site 
and  station  #14,  that  provides  a reasonable  basis  for  correlation 
of  data  from  these  two  locations. 

Similarities  and  differences  between  the  sediment  structure 
at  the  stations  are  apparent  in  short  sections  of  the  reflection 
profiles  (Fig.  3),  and  the  sediment  thickness  in  reflection  travel 
time  is  seen  to  vary  from  about  0.3  to  more  than  2.0  seconds. 

Ridges  and  other  basement  outcrops  in  the  region  are  purposefully 
avoided  (Fig.  3(b)).  Water  depth  for  all  stations  is  about  5100 
meters.  Depth  to  oceanic  basement  obviously  varies  and  at  #8  base- 
ment is  not  detected  below  a deep  and  relatively  strong  reflector. 
Roughness  of  the  sedimentary  horizons  appears  to  be  minimal,  but 
various  degrees  of  roughness  are  exhibited  in  the  basement  returns. 


Seafloor  reflectivity 

Reflectivity  was  measured  in  various  octave  bandwidths 
particularly  at  station  #3  to  see  if  the  determinations  could  be 
refined  (Table  I).  The  average  of  all  reflectivity  coefficient 
measurements  is  0.195,  the  mean  0.209.  There  is  an  upward  trend 
with  increase  in  frequency  which  is  amplified  by  squaring  to  ob- 
tain the  sound  power  reflection  coefficient.  The  reasons  for  this 
trend  are  not  clear  and  cannot  be  resolved  with  the  data  at  hand. 
The  trend  could  quite  possibly  be  caused  by  the  fact  that  the  time 
period  from  the  onset  of  arrivals  which  determines  the  measure- 
ment point  is  held  constant,  while  in  reality  the  effective  pulse 
length  may  become  shorter  with  increase  in  frequency.  Or  the  trend 
may  be  indicative  of  an  effect  of  surficial  reflection  horizon 
spacing  as  a function  of  frequency.  A value  of  0.20  appears  to  be 
a reasonable  representation  of  the  reflectivity  coefficient  for  the 
region. 

A number  of  piston  cores  which  penetrated  the  sediments  1 to 
10  meters  have  been  obtained  in  the  Somali  Basin  environs  during 
four  Woods  Hole  cruiises  (Table  II).  The  upper  sediments  recovered 
were  primarily  calcareous  ooze  or  calcareous  siliceous  ooze  which 
indicates  a degree  of  uniformity  of  the  water-sea  floor  interface 
throughout  the  region.  From  these  samples,  we  anticipate  the  uni- 
formity found  in  the  reflectivity  of  this  interface. 
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Table  II 


Sediment  cores  in  study  area  0°-10°N,  51°-58',E 
(abstracted  from  WHOI  Geological  Samples  Data  File) 


Depth 


Sediment  Description 


Location 


calcareous,  irregular  contacts, 
extensive  mottling, uniform  lithology 
calcareous  ooze,  irregular  contacts 


4499 


calcareous  ooze,  volcanic  glass, 
sharp  contacts,  mottling  uniform 
lithology 

siliceous  clay,  irregular  contacts 


calcareous  ooze,  volcanic  glass, 
uniform  lithology 
calcareous  ooze,  gradational  con 
tacts,  extensive  mottling 
calcareous  siliceous  ooze 


5102 


calcareous  siliceous  ooze,  exten- 
sive mottling 

calcareous  ooze,  gradational  con- 
tacts, cross  bedding 
calcareous  ooze,  gradational  con- 
tacts 

foraminifera 1 sand,  sharp  contacts 

cross  bedding 

calcareous  siliceous  ooze 


4680 


4944 


4864 


calcareous  ooze,  sharp  contacts, 
cross  bedding 

calcareous  siliceous  ooze,  sharp 
contacts,  cross  bedding 
calcareous  ooze,  uniform  lithology 


3697 


4782 


4787 


5104 


calcareous  siliceous  ooze,  CH  100-5 

cross  bedding 

siliceous  ooze,  unfossiliferous  clay.CH  100-4 
cross  bedding 

foraminiferal  sand,  sharp  contacts,  CH  100-4 
cross  bedding 


4576 


We  have  based  some  of  the  later  comparisons  on  this  re- 
flectivity (0.20),  a sound  power  reflection  coefficient  ( (X  ) 
of  0.04,  and  an  effective  pulse  length  of  50  to  80  milliseconds 
with  full  realization  that  these  values  may  vary. 


Reflectance 

Sediment  and  basement  reflectance  in  Table  I are  all  based 
on  an  effective  pulse  length  of  50  milliseconds,  but  on  three  & 's: 
the  standardized  0.040;  the  "average"  0(  at  each  station  (the 
square  of  the  "average  "reflectivity)  ; and  the  measured  in  the 

30-60  Hz  band.  These  are  plotted  in  Figure  12.  Note,  that  these 
reflectance  comparisons  are  made  without  reference  to  impedance 
estimates.  The  choice  of  a 50  millisecond  effective  pulse  length 
as  a basis  for  this  set  of  measurements  causes  the  reflectance  to 
increase  too  rapidly  in  too  many  steps  to  be  accomodated  by  the 
nomograph  in  Figure  7.  A longer  effective  pulse  must  therefore  be 
inferred.  Further,  the  use  of  the  cumulative  energy  curve  pre- 
cludes assessment  of  impedance  inversions — reflected  energy  is 
adaed  on  an  absolute  basis  regardless  of  phase. 

With  these  limitations  in  mind,  this  plot  (Fig. 12)  shows  that 
sediment  reflectance  increases  and  basement  reflectance  decreases 
with  increase  in  sediment  thickness.  The  difference  between  the 
sediment  and  basement  reflectance  measurements  (see  p.21)  for  the  three 
stations  #17,  #14,  and  #15  to  the  south  and  southeast  of  rhe  Basin 
proper  is  consistently  greater  than  at  the  other  three  stations. 

For  shallow  basement  and  thin  sediments  this  can  be  expected.  Yet 
at  #15  the  thickness  of  sediment  is  almost  equal  to  that  at  #4, 
but  the  large  sediment/basement  reflectance  difference  persists  at 
#15  and  does  not  occur  at  #4,  nor  at  #3  and  #8.  This  suggests  that 
there  is  an  entirely  different  depth-distribution  of  impedance  in 
the  southern  part  compared  to  the  northern  part  of  the  region,  a 
difference  that  will  be  more  apparent  when  the  individual  stations 
are  discussed. 

The  general  decrease  in  basement  reflectance  with  increase 
in  sediment  thickness  is  indicative  of  an  impedance  matching 
phenomena  in  which  the  sediments  between  the  low  impedance  water 
and  high  impedance  basalt  act  as  a transmission  coupling  network 
(Appendix  iv) . A chance  distribution  of  impedances  within  the 
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SEDIMENT  THICKNESS  IN  REFLECTION 
TRAVEL  TIME  (seconds) 


Sediment  and  basement  reflectance  as  a function  of  sedi- 
ment thickness  above  oceanic  basement.  (Based  on  effective 
pulse  length  of  50  milliseconds). 

(a)  A general  downward  trend  in  basement  reflectance 
with  increasing  sediment  thickness  is  evident. 

(b)  Sediment  reflectance  does  not  increase  so  rapidly 
with  increase  in  thickness  as  basement  reflectance 
decreases . 

(c)  Difference  between  basement  reflectance  and  sediment 
reflectance  suggests  an  impedance  gradation  in  depth 
at  stations  #4 , #3  and  # 8 that  reduces  basement  re- 
flection and  optimizes  transmission  coupling  into 
basement . 
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sediment  may  satisfy  to  a degree  the  gradation  of  impedance  re- 
quired between  the  two  bodies  of  contrasting  impedance  for 
optimizing  transmission  into  the  basalt  from  the  water.  Apparently 
the  distribution  is  optimized  moreso  at  #4  than  at  #15. 

It  has  been  recognized  (E.T.  Bunce,  G.M.  Purdy,  personal 
communication)  in  deep  sea  refraction  studies  that  the  energy 
from  shots  appears  to  couple  or  transmit  better  into  oceanic  base- 
ment when  sediments  are  present  than  when  not.  But  we  do  not  in- 
tend to  account  for  the  energy  budget  losses  due  to  absorption, 
spherical  spreading,  scattering  etc.  required  to  obtain  a true 
estimate  of  norma]  incidence  transmissivity.  Within  the  framework 
of  our  assumptions  we  can  deal  with  reflectance  and  recognize  its 
relationship  to  transmissivity. 


The  stations 


Figures  13  and  14  show  estimated  impedance  structure  at 
each  station,  the  data  from  which  the  estimates  are  derived,  and 
the  interval  velocity  profile  determined  at  the  location.  The 
section  of  the  normal  incidence  profile  that  the  integrand  repre- 
sents is  indicated;  the  amplitude  squared  trace  is  presented  to 
identify  discrete  reflection  horizons  and  the  basement  return. 
Estimated  impedance  graphs  display  the  maximum  impedance  within 
each  synthetic  layer  computed  from  the  20-percent  energy  steps 
measured  relative  to  the  cumulative  energy  of  the  basement  return. 

In  Figure  13  (stations  #3,  #4,  and  #8)  we  present  impedance 
estimates  based  both  on  percentage  and  on  reflectance.  In  Fig. 

14  the  estimates  based  on  reflectance  are  omitted,  because  they 
are  for  the  most  part  presented  in  Figure  8 (page  24)  and  in 
Figure  15. 


The  most  obvious  differences  in  the  station  data  are  the 
energetic  responses  and  estimates  of  impedances  of  the  sediment 
in  the  two  sets  of  data  as  grouped  in  Figure  13  and  14.  After 
the  initial  return  the  integrands  of  #3,  #4,  and  #8  have  a series 
of  small  steps  up  to  and  including  basement,  whereas  at  #14,  #15 
and  #17  there  is  an  initial  sediment  response  but  little  energy 
cumulated  until  a strong  return  from  basement  is  received.  The 
impedance  steps  are  near  linear  in  the  first  group  and  are  "tread- 
like" in  the  second.  These  differences  seen  on  the  curves  are 
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Stations  #3>  #4,  and  #8  - estimated  impedance  structure  with  depth;  the  seismic 

reflections  from  which  the  estimates  are  derived;  and  interval  velocity. 

(a)  3.5  km  section  of  the  normal  incidence  profile  at  station. 

(b)  Amplitude-squared  trace  of  one  reflection  from  data  window. 

(c)  Average-shape  cumulative  energy  curve  derived  from  reflection  traces  of  data 
window,  normalized  to  a standard  unit  deflection  where  indicated.  Duration 
of  basement  return  indicated  by  bracket. 

(d)  Estimated  acoustic  impedance  based  on  20-percentlle  application  of  model 
(dashed-line)  and  (e)  estimated  impedance  based  on  reflectance  using  a 
50  millisecond  effective  pulse  length  and  a ■ 0.04  (solid  line). 

(f)  Interval  velocity  profile  determined  at  station  (o,  R = critical  refraction 
velocity  determinations). 
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FIG.  14 

Stations  #14,  #15,  and  #17  - estimated  impedance  structure  with  depth;  the  seismic 

reflections  from  which  the  estimates  are  derived;  and  interval  velocity. 

(a)  3*5  km  section  of  the  normal  incidence  profile  at  station. 

(b)  Amplitu ie-squared  trace  of  one  reflection  from  data  window. 

(e)  Average -shape  cumulative  energy  curve  derived  from  reflection  traces  of  data 
window,  normalized  to  a standard  unit  deflection  where  indicated.  Duration 
of  basement  return  indicated  by  bracket. 

(d)  Estimated  impedance  based  on  20-percentile  application  of  model  (dashed-line). 
See  Figures  8 and  15. 

(e)  Interval  velocity  profile  determined  at  station  (o,  R = critical  refraction 
velocity  determinations). 
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much  the  same  as  those  implied  by  reflectance  comparisons  in  Fig. 

12. 


An  intermediate  step  in  the  integrand  occurs  at  stations 
#15  and  #17  between  the  sediment  of  uniform  response  and  the  base- 
ment. Whether  this  transition  should  be  considered  sediment  or 
not  is  unclear.  We  speculate  that  this  may  represent  a transition 
zone  between  the  sediments  and  basalt.  This  structure  is  not  evi- 
dent at  #14.  if  sediment  reflectance  at  #15  is  measured  to  the 
base  of  the  transition,  the  value  for  sediment  reflectance  is  the 
same  as  for  the  base  of  the  sediments  at  #4,  but  the  basement  re- 
flectance is  much  higher  at  #15  than  at  #4  (Table  I).  This  appears 
to  be  a paradox,  but  suggests  that  the  transition  should  be  con- 
sidered either  as  a separate  entity  as  compared  to  the  sediment, 
or  as  a part  of  the  basement  complex. 

In  Figure  15  we  compare  the  impedances  measured  at  DSDP  site 
#235  (Fisher  and  Bunce  et  al,  ibid)  with  the  estimated  impedances 
based  both  on  percentage  and  on  reflectance  at  station  #14.  The 
agreement  between  the  measured  and  estimated  impedance  structures 
is  noteworthy.  For  the  best  fit  in  this  example,  reflectance  is 
based  upon  an  effective  pulse  length  of  80  milliseconds  and  on  CK 
of  .029. 

Station  #14  provides  our  only  comparison  of  the  estimated 
impedance  inferred  from  cumulative  energy  curves  with  the  observed 
structure  obtained  from  samples  at  a specific  site.  Descriptions 
of  surface  cores  at  the  DSDP  site  (Table  II)  obtained  from  R/V 
CHAIN  repeat  findings  near  other  stations,  that  is,  a seafloor  of 
calcareous  or  calcareous  siliceous  ooze.  At  the  drill  hole  cores 
were  recovered  at  eleven  intervals  from  the  seafloor  to  basalt  650 
meters  below.  The  core  at  580  meters  encountered  moderate  brown 
clay  and  recovered  only  3%  of  the  material  drilled,  substantially 
less  than  the  other  intervals.  This  change  in  lithology  and  com- 
petence may  be  indicative  of  the  transition  between  sediment  and 
basement  seen  in  the  integrands  at  stations  #15  and  #17,  but  such 
a transition  is  not  apparent  at  #14. 

The  report  (ibid)  gives  the  acoustic  impedance  derived  from 
velocity  and  density  measurements.  These  are  not  in  situ  velocity 
measurements  so  a value  such  as  1.47  km/sec  for  nanno  ooze  near 
the  seafloor  may  be  less  than  if  measurements  were  made  at  the 
temperature  and  pressure  at  depths. 
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Our  impedance  estimates  versus  travel  time  together  with 
the  impedance  measurements  from  DSDP  site  #235,  appear  to  con- 
firm the  intuitive  interpretation  of  the  widely  observed  "trans- 
parent layer"  phenomenon.  The  agreement  is  encouraoing  in  both 
the  percentage  and  reflectance  applications  of  the  lov/- frequency 
reflection  model. 


Interval  velocity  structure  re  impedance  structure 

Of  the  two  groups  of  stations,  those  to  the  south  (Fig.  14) 
show  the  more  orderly  relationship  between  impedance  structure  and 
interval  velocity  structure.  Both  velocity  and  impedance  increase 
together  in  a logical  fashion.  Impedance  estimates  based  on  re- 
flectance and  the  sound  power  reflection  coefficient  also  run  close 
to  those  based  on  percentage  (Figs.  8 and  15) . 

In  the  northern  group  (Fig.  13)  we  have  included  estimates 
of  impedance  based  on  reflectance  and  the  sound  power  reflection 
coefficient  (50  ms  effective  pulse  length  and  CK , = .040)  as  well 
as  those  based  on  percentage.  Impedance  curves  based  on  re- 
flectance and  percentage  run  close  to  each  other  onlv  at  station 
#4. 


At  stations  #3  and  #8  there  are  pronounced  velocity  inversions, 
which  conceivably  would  cause  stronger  impedance  inversions,  than 
those  at  #4.  In  our  model  and  in  our  method  of  determining  re- 
flectance we  assume  no  inversions.  Squaring  removes  phase  infor- 
mation thus  making  all  contributions  to  reflectance  positive  con- 
tributions. Further  when  using  the  percentage  approach  the  impedance 
estimates,  which  are  maximum  values,  are  made  to  fit  between  the 
impedance  of  seawater  at  depth  and  that  of  basalt.  Therefore,  it 
can  be  seen  that  the  impedance  estimates  based  on  reflectance  would 
increase  with  depth  more  rapidly  than  both  the  actual  impedance 
and  the  percentage  based  estimates  if  there  were  appreciable  in- 
versions. Evaluation  of  the  impedance  estimate  must  therefore  be 
partially  based  on  a knowledge  of  the  velocity  structure.  Con- 
versely when  the  reflectance  based  impedance  estimate  increases  too 
rapidly,  the  presence  of  velocity  inversions  such  as  we  measured 
in  the  sonobuoy  data  is  suspected. 
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Graphic  comparisons  of  reflectance  at  the  stations 

Expanded  scale  integrands  of  the  six  stations  are  com- 
pared in  terms  of  reflectance  and  a common  sound  power  reflection 
coefficient  in  Figure  16.  Although  these  comparisons  can  only 
be  made  on  an  approximate  basis  because  of  the  uncertainties  in 
pulse  length  and  d , the  comparison  brings  out  similarities  and 
differences  more  clearly  than  in  Figures  13  and  14.  The  remark- 
able similarity  in  the  sediment  response  at  station  #3  and  #4  cor- 
roborates what  might  be  expected  at  adjoining  stations.  At  station 
#8,  compared  to  #3  or  #4,  a more  gradual  rise  at  the  very  beginning 
of  the  integrand  after  the  onset  suggests  a change  in  seafloor 
character  that  accounts  for  the  apparent  decrease  in  seafloor  re- 
flectivity at  #8  (Table  I).  The  similarity  of  the  acoustic  strati- 
graphy in  the  first  1.25  sec.  below  seafloor  seen  on  the  reflection 
profiles  of  stations  #3  and  #8  (Figure  3a)  is  less  evident  in  the 
energetxc  response.  There  are  a few  events  during  the  sediment 
returns  that  appear  to  correlate,  but  the  overall  rate  of  increase 
of  the  integrand  at  #8  is  considerably  less  than  at  #3  and  much 
closer  to  that  at  #14  and  #15.  This  suggests  that  the  transition 
in  impedance  structure  between  #3  and  the  southern  group  of 
stations  (Fig.  14)  begins  to  occur  at  #8  only  150  km  southwest  of 
#3. 


Contoured  traverse  profiles 

Although  our  short  sections  of  contoured  data  show  much  the 
same  reflectance  structure  as  is  found  at  the  respective  stations, 
it  is  of  interest  to  compare  the  different  contouring  methods  and 
to  speculate  as  to  continuity  or  lack  of  continuity  of  seabed 
characteristics  along  the  profiles. 

In  Figures  9 and  10  (pgs.  26,  27)  successive  integrands  and 
contours  of  the  integrands  along  a traverse  are  presented  to 
examine  the  continuity  of  reflectance  along  the  traverse  as  well 
as  to  demonstrate  the  analysis  process.  The  waveforms  shown  in 
Figure  10  are  taken  from  the  very  beginning  of  the  profile  at 
station  #3  shown  in  Figure  17.  Data  in  Figure  17  are  not  smoothed 
to  the  same  extent  by  the  running  average  algorithm  demonstrated 
in  Figure  10,  but  are  otherwise  based  on  the  same  analysis  pro- 
cedures. 
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REFLECTION  TRAVEL  TIME 
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Comparison  of  seabed  reflectance  at  the  six  stations 
using  the  average-shape  integrand  from  each  station, 
normalized  to  a standard  unit  deflection  at  50  milli 
seconds  after  onset  of  returns.  Basement  returns 
indicated  by  brackets. 


reflection  travel  time  in  seconds 


STATION  *3 


STATION  *14 


Seismic  reflection  profiles  at  station  #3  and  #14  with  overlays  of  contoured 

cumulative  energy. 

(a)  #3  profile  crosses  a thick  (2.3  second)  section  of  sediments;  basaltic 
basement  is  the  rough  faint  horizon  occurring  between  8.9  and  9.2  seconds 
travel  time. 

(b)  #14  profile  crosses  a thin  (0.7  second)  section  of  sediments;  basaltic 
basement  is  the  rough  strong  horizon  occurring  at  7.2  to  7.6  seconds  travel 
time . 

(c)  The  reflections  occurring  between  the  seafloor  reflection  at  6.8  seconds 

and  the  interface  at  7.6  seconds  in  profile  "a"  appear  to  be  markedly  uniform 
as  in  Fig.  10.  Normalizing  each  succeeding  cumulative  energy  curve  to  a 
common  deflection  at  7.6  seconds  travel  time,  the  succeeding  curves  are  then 
contoured. to  show  the  travel  times  of  the  points  on  the  cumulative  energy 
curve  having  0.2,  0.4,  0.6,  0.8,  and  1.2  of  the  unit  deflection  at  7.6 
seconds.  Assuming  a constant  impedance  to  exist  at  the  interface  at  7.5  to 
7.6  seconds  travel  time,  these  contours  serve  to  show  the  lateral  fluctuations 
in  the  impedance  of  the  overlying  sediment  along  the  traverse. 

(d)  The  cumulative  energy  curves  can  be  examined  as  in  Fig.  10  by  assuming  a 
uniform  reflectivity  of  the  seafloor  Interface  along  the  traverse  and  then 
sectioning  the  deflection  of  the  cumulative  energy  curve  to  determine  the 
travel  times  at  which  Integer  multiplies  of  the  initial  deflection  occur. 
Irregularities  in  travel  time  Indicate  there  is  fluctuation  in  the  initial 
seafloor  return  even  in  areas  of  uniform  sedimentation. 

(e)  The  procedure  outlined  in  "c"  is  applied  to  profile  "b"  using  oceanic  base- 
ment as  the  reference  horizon  and  using  a norm  point  at  7.6  seconds,  about 
0.1  second  after  the  start  of  the  basement  return. 

(f)  Using  the  same  procedure  in  "d",  the  sediment  sequence  is  shown  to  be  of 
relatively  uniform  reflectance  and  hence,  Impedance.  This  Inference  is 
clearer  here  than  in  profile  "e"  since  this  contouring  procedure  is  not 
affected  by  the  topography  in  the  basaltic  basement. 


In  Figure  17a  the  evident  uniformity  of  the  upper  sediments 
seen  in  the  conventional  reflection  profile  across  station  #3  is 
substantiated  by  the  contoured  reflection  energetics  (Fig.  17c) . 

The  uniformity  is  best  displayed  using  the  consistent  reflection 
horizon  at  depth,  7.6  seconds  travel  time  (about  0.8  secs,  below 
the  seafloor)  as  the  norm  point  for  contouring.  Contours  based 
upon  integer  multiples  of  an  assumed  constant  sound  power  re- 
flection coefficient  fluctuate  more,  but  show  relatively  consistent 
or  coherent  shifts  across  the  contours  'Fig.  17d).  This  indicates 
that  seafloor  reflectivity  here  varies  from  shot-to-shot  and  varies 
more  than  the  percentages  of  reflectance  from  the  seafloor  down 
to  the  uniform  reflector  at  depth.  There  is  however  a uniformity 
in  the  reflectance  contours  (Fig.  17d)  that  is  evident  in  spite  of 
the  oscillations. 

In  contrast,  at  station  #14  (Fig.  17e)  the  40%  contour  shows 
considerable  variation.  Small  changes  in  the  reference  energy 
from  basement  causes  this  contour  to  meander  greatly  because  the 
reflectance  and  impedance  of  the  sediments  are  so  uniform.  Re- 
flectance contours  that  are  related  to  seafloor  reflectivity  at 
this  station  (Fig.  17f)  indicate  better  the  continuity  of 
structure  and  suggest  a more  uniform  seafloor  reflectivity  than 
at  station  #3.  (A  continuous  profile  that  crossed  from  a region 
as  at  #3  into  a region  as  at  #14  would  have  been  most  interesting! ) 


Mapping  characteristics 

Too  few  measurement  stations  are  included  in  this  study  to 
warrant  mapping  different  characteristics.  Changes  in  the  dif- 
ferent parameters,  reflectivity,  reflectance,  time  constants,  and 
Area  would  be  evident  in  traverses  from  one  physiographic  region 
to  another,  but  our  data  here  are  limited  to  the  relatively  short 
tape  recordings  taken  during  the  sonobuoy  stations. 


CONCLUSIONS  & RECOMMENDATIONS 

1.  Our  major  conclusion  is  that  an  "audit”  of  the  seabed  re- 
flected energy  (reflectance)  yields  impedance  estimates  commensu- 
rate with  the  available  drilling  samples,  the  compressional  wave 
velocities  determined  from  wide-angle  reflection-refraction  measure- 
ments, the  descriptive  appearance  of  the  normal-incidence  profiling 
recordings  and  the  seafloor  sampling. 
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2.  Two  distinct  impedance-depth  structures  are  evident  in 

the  subject  environs  of  the  Somali  Basin:  (a)  A near-linear 

increase  with  depth  at  stations  #3,  #4  and  #8.  (b)  A non-linear 

tread-like  increase,  having  a relatively  uniform  low  impedance 
throughout  the  sediment  column  down  to  a transition  zone  or  to 
basement  at  stations  #14,  #15  and  #17.  Estimated  impedance  pro- 
files based  on  seafloor  reflectance  show  the  same  trends.  The 
discrepancy  between  estimates  based  on  percentage  and  reflectance 
at  stations  #3  and  #8  can  be  explained  by  accounting  for  the 
presence  of  impedance  inversions  by  the  pronounced  interval  velo- 
city inversions  measured  at  the  stations.  Otherwise  differences 
are  probably  the  result  of  uncertainties  in  seafloor  reflectivity 
and  pulse  length.  Agreement  between  impedance  estimates  based 

on  percentage  and  reflectance  is  quite  acceptable  when  these  two 
parameters  are  adjusted  to  the  pulse  length  and  sound  power  re- 
flection coefficient  measured  at  30  to  60  Hz. 

3.  Seafloor  reflectivity  in  the  30-60  Hz  passband  is  quite 
uniform  at  the  six  stations,  median  0.19,  but  the  subtle  dif- 
ferences observed  appear  to  be  the  key  to  agreement  between  the 
percentage  and  reflectance  based  estimate  of  impedance.  When 
examined  over  a wider  range  of  frequencies,  the  reflectivity  shows 
more  variability,  from  0.12  to  0.30.  Seafloor  cores  consist  pri- 
marily of  calcareous-siliceous  ooze,  suggesting  uniformity  in  sea- 
floor composition  as  well  as  reflectivity. 

4.  Since  reasonable  agreement  is  obtained  between  the  cored  re- 
sults at  DSDP  #235  and  the  estimates  of  acoustic  impedance,  it 
appears  that  the  assumptions  r£  internal  reflection  multiples, 
plane  wave  propagation,  absorption  and  scattering  are  practical 
assumptions  to  make  when  working  at  low  frequencies  such  as  the 

30  to  60  Hz  passband. 

5.  We  conclude  that  there  is  need  for  further  application  of 
this  type  of  analysis  and  that  measurements  of  reflectivity,  re- 
flectance and  impedance  estimates  are  of  unique  utility  in  inter- 
pretation of  seismic  reflection  profiles.  The  application  of  the 
analysis  will  be  of  specific  interest: 

a.  to  provide  better  assessment  of  the  characteristics  and 
the  continuity  of  the  characteristics  of  the  seabed  for  implanting 
materials,  instruments  and  drilling.  This  type  of  assessment  is 
germane  to  better  correlation  between  the  acoustic  record  and  the 
characteristics  of  cored  sediments  and  to  an  improved  evaluation 
of  the  lateral  homogeneity  of  materials  between  core  or  drill  holes 
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b.  to  examine  within  the  transitional  phenomena  along  pro- 
files the  manner  in  which  seafloor  reflectivity  relates  to  changes 
in  grain  size  and  sorting,  and  in  differentiation  between  pelagic, 
turbidity,  and  contourite  deposits.  This  requires  more  resolution 
than  available  in  this  study  but  could  be  accomplished  by  using 
higher  frequencies,  i.e.  sound  pulses  with  center  frequencies  at 
3.5  and  12  kHz.  it  is  relatively  straight  forward  to  include 
estimates  of  transmission,  absorption  and  scattering  losses  at 
the  higher  frequencies- 

c.  to  provide  additional  information  with  which  to  evaluate 
determinations  of  interval  velocities  from  wide  angle  reflection- 
refraction  profiles,  from  sonobuoys  or  multi-channel  arrays.  The 
reflectivity  of  the  seafloor  restricts  the  range  of  velocities  of 
the  first  layer,  and  the  reflectance,  the  impedance  of  the  deepest 
observed  reflection.  Should  the  impedance  estimates  increase  to 
an  impedance  value  such  as  that  of  basalt  at  a travel  time  earlier 
than  observed,  the  presence  of  velocity  and/or  impedance  inversions 
is  suggested. 

d.  to  make  quantitative  measurements  of  the  shot-by-shot  re- 
flection variability,  to  examine  and  possibly  to  identify  the 
mechanisms  producing  this  variability,  and  to  use  the  variability 
as  a measure  of  seafloor  and  seabed  character. 
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Appendix  I Reflectance  of  a seabed  consisting  of  two  layers. 

A simplified  representation  of  ray  paths  for  seismic 
profiling  is  shown  in  Figure  1-1.  An  intermediate  layer  of 
material  with  a characteristic  impedance,  P,  C, , is  between  the 
water  (fpC-p)  and  the  acoustic  basement  material  ( ) . a 

continuous  acoustic  wave  with  an  initial  amplitude,  ^ strikes 
the  first  interface  at  normal  incidence.  (For  clarity  the 
figure  is  drawn  with  the  rays  at  an  angle.)  At  the  first 
interface  part  of  the  beam  is  transmitted  and  part  reflected. 
The  transmitted  beam  encounters  the  second  interface  where 
it  is  partly  reflected  and  partly  transmitted.  The  reflected 
beam  returns  to  the  first  interface  where  again  reflection 
and  transmission  occur.  This  process  continues  with  an  in- 
finite number  of  reflections  and  a steady-state  situation  is 
obtained . 

At  the  first  interface  the  reflection  coefficient  from 
water  to  the  intermediate  layer  is 


Z,  - / 

i,  + / Eq.  1-1 

where  is  the  ratio  of  the  acoustic  impedances,  /°<  c' //4?o  . 
For  a reflection  at  this  interface  within  the  layer,  the  re- 
flection coefficient  r1Q  is  equal  to  ( -rQ:L)  . The  trans- 
mission coefficient  across  the  first  interface  from  the  water 
to  the  intermediate  layer  is 
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Eq.  1-2 


But  the  transmission  coefficient  across  the  first  interface 
from  the  layer  to  the  water  is 


1 interface 
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The  power  reflection  coefficient,,^,  is  the  ratio  of  the 
intensity  of  the  reflected  wave  to  the  intensity  of  the  in- 
cident wave.  And  since  both  are  in  the  same  medium,  the 
power  reflection  coefficient  becomes  the  ratio  of  the  squares 
of  the  amplitudes  or 
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Eq.  1-4 


The  power  transmission  coefficient,  T,  is  the  ratio  of 
the  intensity  of  the  transmitted  wave  to  the  intensity  of  the 
incident  wave.  However  in  this  situation  the  two  waves  are 
in  different  media.  Thus  the  power  transmission  coefficient 
is 
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Eq.  1-5 
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Notice  that  T^  is  equal  to  the  product  of  tQl,  and  t1Q 
that  Of,  + Tx  = 1. 


and 


In  the  more  common  situation  the  intensity  of  the  trans- 
mitted waves  is  computed,  but  for  seismic  profiling  the  ob- 
served quantity  is  the  reflected  waves.  Although  the 
instantaneous  amplitude  is  recorded,  the  energy  (intensity) 
is  the  quantity  that  is  used  in  the  analysis.  According  to 
the  Principle  of  Superposition,  the  amplitudes  of  the  re- 
flected components  must  be  added  with  their  proper  phase. 
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Then  the  intensity  is  computed  from  the  resultant  amplitude. 

The  phase  change  due  to  reflection  at  an  interface  is 
assumed  to  be  0 or  180  degrees.  The  phase  change  t t due 
to  the  path  difference  is  (path  difference/wave  length)  or 
^ where  h is  the  thickness  of  the  layer.  (Note  that 
in  computing  the  intensities,  phase  changes  of  180°  will  not 
effect  the  results.) 

In  Figure  1-1  the  amplitudes  of  the  waves  are  given  after 
each  reflection  or  transmission,  thus  amplitude  of  the 

first  reflected  wave  B B.rg^  and  the  first  transmitted 

wave,  Bq t q ^ . The  order  jf  the  subscripts  is  related  to  the 
direction  of  wave  at  each  event.  The  phase  changes  are  not 
included  in  this  figure. 

The  component  that  enters  the  layer  is  reflected  from 
the  second  interface,  and  finally  it  re-enters  the  water. 

So  its  amplitude  is  BQtQir^tlo  * The  resultant  amplitude 
including  the  phase  difference  between  successive  waves  can 
be  written  as 

BrW=  B.L0|v;2-t,oei')-v  Y>'Lo,r'z  ejl  H . .. 


The  common  terms  may  be  factored  from  the  infinite  series  so 


But  = T^ . Also  the  series  in  the  brackets  is  a geo- 
metric progression  with  a ratio  (r^r.-e  ),  and  thus  it  is 

readily  summed.  The  reflected  amplitude  may  thus  be  written 
as 


This  equation  may  be  simplified  if  all  quantities  are  over 
the  common  denominator.  In  the  algebraic  rearrangement,  the 
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identities  = 1 - and  rQ1  =-r1Q  are  used.  The  equation 

may  then  be  written  as 


I • t'x  l_ 

\ ' r \r  \ 

12-10  _J 

The  resultant  amplitude  is  a complex  quantity.  Thus 
the  intensity  is  the  product  of  the  resultant  amplitude  and 
its  complex  conjugate.  By  using  the  Euler  Relation  the  re- 
flected intensity  becomes 
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Or  in  terms  of  the  power  reflection  coefficients 


1 


llg 

I 


Irefl/I| 


0 


o< , / <*z  — x 


/ c>( , [X  ^ 2 ^1 fp(  t M i CdS  c/" 


Eq.  1-6 


which  is  what  we  have  termed  the  reflectance,^*  , for  the 
single  layered  seabed  with  two  interfaces. 
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Appendix  II  Comparison  of  the  reflectance  computations 
for  a two  layer  seabed  as  a function  of  ( ’ 


The  equation  for  the  reflectance  derived  in  Appendix 
I is  for  a steady  state  wave  with  a single  frequency.  The 
signals  used  in  seismic  profiling  are  broadband  pulses. 

The  electronic  processing  of  the  received  signal  permits 
the  analysis  to  be  done  in  a smaller  band  of  frequencies 
but  not  on  a single  frequency.  Also  the  computational  prob- 
lems involved  in  using  Eq.  1-6  can  be  reduced  by  the  use  of 
an  approximation  for  a broader  band  signal.  The  reduced  form 
for  the  reflectance  from  a single  layer  with  two  interfaces 
is  given  by  Eq.  3. 


= at.  + T«z 


Eq.  II-l 


l 

H 


The  range  of  acoustic  impedances  commonly  encountered 
in  deep  sea  seismic  profiling  is  from  water  to  basalt  or  a 
ratio  with  respect  to  water  of  about  1 to  10.  By  picking 
values  of pC  between  water  and  basalt  the  reflectance  can  be 
computed  by  Eq.  II-l  and  by  Eq.  1-6.  However  in  the  case  of 
Eq.  1-6  the  effective  path  difference  must  be  included.  So 
for  the  response  due  to  a relatively  broadband  signal  we  will 
compute  the  reflectance  by  Eq.  1-6  for  path  differences  be- 
tween 0°  and  180°  in  steps  of  20°.  Then  we  will  take  the 
average  for  comparison  purposes. 

Table  II-l  gives  the  results  of  these  computations. 

Here  Z ^ /fipo  is  the  relative  acoustic  impedance  between 

the  layer  and  water.  By  letting  Z^  vary  between  1 and  10, 
the  effect  is  to  consider  the  layer  between  the  water  and 
the  basalt  to  vary  in  acoustic  impedance  from  water  to  basalt. 


Table  II-l 


Reflectance 


Z1 

Eq.  II-l 

Eq.  1-6 

1 

.67 

.67 

2 

.46 

.47 

4 

.44 

.43 

6 

.53 

. 52 

8 

.61 

.60 

10 

.67 

.67 

small 

The  differences  shown  in  this  table 
, certainly  within  the  experimental 

are  relatively 
error  entailed  in 

the  measurements. 
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Appendix  III 


A simplified  derivation  for  the  reflectance 
from  a many  layered  seabed. 


In  Appendix  I the  reflectance  from  a seabed  with  one 
layer  and  two  interfaces  is  analyzed.  This  analysis  sums 
the  amplitudes  of  all  of  the  multiple  reflections  with  the 
proper  phase  included.  The  relative  intensity  is  then 
obtained  from  the  square  of  the  relative  amplitude. 


In  Appendix  II  we  show  that  the  relative  intensity  as 
computed  from  the  average  over  the  acoustic  wavelengths  for 
a broad  band  signal  is  approximately  the  same  as  the  relative 
intensity  obtained  by  neglecting  the  phase  relations  between 
successive  reflections.  In  the  situation  most  commonly  en- 
countered in  the  seabed,  the  successive  layers  tend  to  be 
compressed  with  depth  producing  increasing  values  of  acoustic 
impedance.  In  other  words,  the  power  reflection  coefficient 
between  the  individual  layers  is  small  and  the  transmissivity 
is  large. 

These  two  conditions,  averaging  of  the  phases  and  the 
small  power  reflection  coefficient  between  layers,  permits 
a simplified  derivation  for  the  reflectance  for  a multi- 
layered seabed.  Since  we  are  not  considering  the  phase  re- 
lations, we  may  deal  directly  with  the  intensities.  And 
since  the  power  reflection  coefficient  is  small  we  need  to 
consider  only  the  first  reflection  from  each  surface.  With 
these  conditions,  the  reflection  process  is  indeed  that  shown 
in  Figure  III-l.  (Note  that  the  incoming  and  reflected 
signals  are  at  normal  incidence  but  they  are  drawn  at  an  angle 
to  clarify  the  process) . The  intensity  of  the  incoming  signal 
is  iQ.  At  the  first  interface,  it  is  partly  reflected  (c?<,  Iq) 
and  partly  transmitted  (T^Iq).  When  the  transmitted  signal 
encounters  the  second  interface,  partial  reflection  (T^p(?I  ) 
and  partial  transmission  occurs*  The  upgoing  signal 

returns  to  the  first  interface  where  the  transmitted  portion 
(Tf  (X ' Iq)  is  the  part  that  returns  to  the  receivers  on  the 
surface.  The  process  continues  at  the  successive  surfaces 
as  shown  in  the  diagram. 


To  show  that  the  second  reflection  term  is  much  smaller 


warn 


T1  T2T3T4I0 


A schematic  diagram  to  show  the  reflection  of  a 
pulse  from  a layered  seabed.  (The  incoming  pulse 
is  at  normal  incidence.  The  figure  is  drawn  at 
an  angle  to  clarify  the  reflection  paths). 
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and  that  it  may  be  negligible  for  the  accuracy  of  our  measure- 
ments, we  will  compare  the  magnitude  of  the  first  and  second 
reflections  from  the  first  interface.  The  intensity  of  the 
first  reflection  is,  as  described  above,  (T-^2  p(z  Iq).  For  the 
signal  that  is  reflected  the  second  time  from  the  second  inter- 
face before  returning  to  the  receiver,  the  intensity  is  ^ ^ q) 
or  ( o(to(  ) times  smaller  than  the  first  reflection.  With  ot, 
and  p(z  on  the  order  of  one-tenth,  their  product  is  on  the  order 
of  one-hundredth  or  negligibly  small. 

Thus  the  signal  returned  to  the  receiver  from  a specific 
interface  is  the  product  of  the  square  of  the  power  transmission 
coefficients  for  the  intervening  interfaces  and  the  power  reflection 
coefficient  for  that  interface.  These  return  signals  are  shown 
for  a four  layered  seabed  with  five  interfaces  in  Figure  III-l. 

Since  the  reflectance  for  the  seafloor  is  the  ratio  of  the  net 
reflected  signal  to  the  incident  signal,  the  reflectance  is  readily 
calculated  as: 


rXr;^ 


Eq.  III-l 


which  is  the  same  as  Eq.  4. 

In  Figure  III-2  the  reflectance  for  the  received  signal  for 
three  types  of  seabeds  is  plotted.  These  reflectances  are  com- 
puted from  Eq.  III-l  and  are  mathematical  models.  They  are  not 
intended  as  an  attempt  to  model  data?  rather  these  figures  are 
intended  to  show  the  effects  of  different  acoustic  impedances  for 
the  successive  layers.  The  first  curve  is  the  result  of  equal  in- 
crements in  the  relative  acoustic  impedances  between  water  and 
basalt.  In  the  second  curve,  the  effect  of  a strong  reflector  at 
the  second  interface  is  shown.  For  the  third  curve  the  interfaces 
reflect  approximately  equal  energies. 


WHEN  Z 


Computed  energy  integrands  for  selected  relative  impedances  of 
the  seafloor. 


Equal  increments  in  the  relative  acoustic  impedances  between 
water  and  basalt  Z„  = 1,  2,  4.  6,  8,  10. 


The  effect  of  a strong  reflector  at  the  second  interface 


For  interfaces  that  reflect  approximately  equal  energies 
Z„  = l,  1.5,  2.4,  3-7,  6.0,  10. 
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Appendix  IV.  Gradation  of  seabed  acoustic  impedance  with 
depth  re  reflectance. 


The  relationship  of  the  magnitude  of  reflectance  at  various 
travel  times  on  the  integrand  to  the  magnitude  of  the  return  from 
a reference  horizon  such  as  basaltic  basement  gives  an  indication 
of  the  impedance  characteristics  with  depth.  Figure  IV-1  indi- 
cates the  reflectance  from  a sequence  of  interfaces  between 
synthetic  layers,  each  of  which  produces  an  equal  reflected  energy 
increment.  The  maximum  reflectance  is  constrained  by  an  upper 
limit  of  impedance,  that  of  basalt.  If  the  seafloor  were  basalt, 
the  Rayleigh  reflection  coefficient  would  be  approximately  0.83 
and  about  69%  of  the  incident  energy  would  be  re|lecte^  back  into 
the  water,  using  acougtic  impedances  of  1.6  x 10  kg/nr  sec  for 
sea  water  and  18  x 10  for  basalt. 


This  is  apparent  in  Figure  IV-1  when  entering  the  graph  at 
the  appropriate  reflectivity.  One  interface  is  involved  and  it 
is  characterized  by  a high  impedance  contrast.  As  more  inter- 
vening layers  are  introduced  between  the  seawater  and  basalt,  each 
having  appropriate  impedances  to  fit  the  equal  energy  step  model, 
the  reflectance  in  relation  to  basaltic  basement  decreases.  Neg- 
lecting transmission  losses,  absorption  and  scattering,  the  total 
amount  of  energy  reflected  back  into  the  water  is  reduced  and  the 
amount  of  energy  transmitted  into  the  basalt  presumably  increases. 


This  same  effect  may  also  be  examined  by  plotting  the  re- 
flectance of  a three-layer  geological  model  (sea-water,  inter- 
mediate-layer, basalt)  in  which  the  impedance  of  the  intermediate 
layer  varies  from  that  of  seawater  to  that  of  basalt  (Fia.  IV-2) . 


Carrying  this  examination  one  step  further  to  a four  layer 
model  bounded  by  seawater  and  basalt  with  two  intermediate  layers, 
an  approximate  plot  of  "total"  reflectance  contours  computed  re- 
lative to  the  basaltic  basement  return  shows  the  effects  of  an 
impedance  inversion  (Fig.  IV-3).  if  the  impedance  of  the  third 
layer  is  low  relative  to  that  of  basalt  as  well  as  that  of  the 
overlying  layer  (an  impedance  inversion  within  the  intermediate 


energy 


reflectance 

(%  of  incident  energy) 


.1  .2  .3 

REFLECTIVITY  OF  INTERMEDIATE  LAYER 


FIG.  IV-2  Plot  showing  change  in  reflectance  as  a function  of 

the  impedance  of  one  intermediate  layer.  Total  amount 
of  energy  reflected  back  into  the  water  is  reduced 
when  we  introduce  a layer  of  intermediate  impedance 
between  the  water  and  basalt. 


IMPEDANCE  OF  SECOND  LAYER 
x 106  ray/s  ( kg/m2sec ) 
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15  3.0  4.5  6.0  7.5  9.0  10.5  12.0  13.5  15.0  16  5 

IMPEDANCE  OF  THIRD  LAYER 
xIO6  ray/s  (kg/m2sec) 


FIG.  IV-3  Plot  showing  change  In  reflectance  as  a function  of  the 
Impedance  of  two  Intermediate  layers  between  water  and 
basalt . 


layers)  the  reflectance  is  seen  to  increase. 

These  examples  in  conjunction  with  the  analysis  of  data 
that  are  presented  illustrate  ways  in  which  an  improved  under 
standing  of  seismic  reflection  energetics  can  lead  to  refined 
qeoioqical  interpretations  of  the  seismic  reflection  profile. 


A 


1 
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Appendix  v.  An  Estimate  of  the  Effect  of  Noise  on  the 
Cumulative  Energy  Curves. 

Two  "signals"  are  arriving  simultaneously  for  analysis, 
the  desired  signal  is  a coherent  time  series  S(t)  while  the 
noise  is  a random  time  series  N(t).  For  the  cumulative 
energy  curve  we  want  the  sum  of  the  energy  from  the  start  of 
the  signal  up  to  a time  T,  or  JJ'fsU^dt  . since  both  signal 
and  noise  are  recorded,  their  instantaneous  amplitudes  add. 

So  the  signal  to  be  processed  is 


The  first  term  is  the  desired  cumulative  energy  of  the 
acoustic  signal.  The  second  term  in  the  integral  is  the 
cumulative  energy  of  the  noise  in  the  interval.  We  postu- 
late that  the  received  noise  background  is  sufficiently  uni- 
form that  we  can  use  an  average  noise  signal.  The  average 
noise  level  is  computed  in  a 500  millisecond  window  before 
the  first  seafloor  return,  and  may  be  designated  as  ((*>(£ Jj  . 
Thus  the  second  term  may  be  expressed  as  ))*</{  -((nMY)T • 

Actually  in  the  digital  processing  of  the  data  the  average 
noise  energy  per  digital  interval  is  subtracted  from  the 
total  energy. 

The  third  term,  Ja  2 s(t)  N(t)  cLt  , is  more  difficult  to 
assess.  It  is  the  integral  of  the  product  of  two  uncorrelated 
time  series  so  the  instantaneous  magnitude  of  this  term  is 
determined  by  the  phase  relations  between  S(t)  and  N(t).  Both 
S ( t ) and  N(t)  have  positive  and  negative  values  and  their 
product  may  also  have  both  positive  and  negative  values.  if 
this  term  is  summed  over  a long  time  its  value  should  be  zero, 
however,  the  50  to  100  millisecond  pulse  duration  of  the 
seismic  source,  the  expected  length  of  a reflection  event, 
is  not  long  enough. 

The  result  of  the  digital  processing  of  the  data  should 
thus  produce  the  following: 

X f ‘(sil  » V*  + xjo  5(i) 


Obviously  the  lower  bound  of  the  cross-product  term  is 
zero.  But  zero  is  not  its  expected  value  in  a short  time 
interval,  rather  this  term  is  indeterminabxe  although  small. 

For  the  first  seafloor  arrival  with  a signal-to-noise 
ratio  greater  than  20  dB,  the  term  JT(s(t\)lcL'k  will  be  so 
much  larger  than  the  cross-product  term  that  its  effect  is 
trivial.  But  for  the  second  arrival  with  a signal-to-noise 
ratio  near  0 dB,  the  magnitudes  of  S(t)  and  N(t)  are  nearly 
equal.  In  this  case  the  cross-product  term  is  very  likely 
to  be  sufficiently  large  so  that  it  cannot  be  neglected  in 
comparison  with  cLt  . Thus  it  introduces  a random 

error  in  the  cumulative  energy  of  the  second  arrival  and 
also  in  the  computation  of  either  the  sound  power  reflection 
coefficient  ( cxf  ) or  the  Rayleigh  reflection  coefficient  (r)  . 
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